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Preface 


This report contains the design of both a Dameus and a 
Giromill for use on Mars. The report has been organized so that the 
interested reader may read only about one machine without having 
to read the entire report. Where components for the two machines 
differ greatly, separate sections have been allotted for each machine. 
Each section is complete; therefore, no relevant information is 
missed by reading only the section for the machine of interest. Also, 
when components for both machines are similar, both machines have 
been combined into one section. This is done so that the reader 
interested in both machines need not read the same information 



i. 

L 

[ 


Table of Contents 


Item 

Statement of Purpose 
Environmental Concerns 
Wind Speed Design Decisions 
Wind Turbine Selection 
Darrieus Blade Design 
Introduction 

Size and Tip-Speed Calculations 
Troposkien Shape and Optimization 
Airfoil Selection 
Material Selection 
Structural Analysis; Mass of Blades 

Giromill Blade Design 

Dimensions and Tip-Speed 
Strut Positioning 
Blade Structure 
Connections to Struts 
Stresses and Factor of Safety 
Total Mass 
Airfoil Selection 
Material Selection 

Shaft Design 

Gear Design, Lubrication of Gears and Bearings 
Generator and Electronics 
Deployment 
Conclusions 
Future Work 
Machine Configurations 
Darrieus 
Giromill 

Deployment Illustrations 

Accurate Troposkien Shape 

Illustration of Generator and Electronics 

References 

Appendices 


Page 

1 

2 

3 - 4 
5 - 6 

7 

8 

9-14 
15 - 18 
19 

20 - 24 

25 - 26 

26 - 27 

27 - 28 

28 - 29 

29 - 30 
30 

3 1-34 
35 
36 - 4 1 
42 
43-45 
46-47 

48 

49 

50 
5 1 

52 - 55 

56 

57 
58 - 59 

see Appendix 
Table of Contents 



List of figures 


Figure ^*5 

Figure 0 : Troposkien Notation T10 

Figure 1 : Free-Body Diagram of Troposkien T 1 1 

Figure 2 : Power-to-Mass for Darrieus 

in main report 1 1 

in appendix (larger version) T37 

Figure 3 : Shape Comparisons for Darrieus 

in main report 1 2 

in appendix (larger version) T38 

Figure 4 : Tension and Stress Ratios 

in main report 1 4 

in appendix (larger version) T39 

Figure 5 : Coefficient of Power 

for Darrieus 1 6 

for Giromill 3 3 

Figure 6 : Incident Velocity vs. Free Stream Velocity 

for Darrieus 1 2 

for Giromill 3 4 

Figure 7 : Coefficient of Moment 

for Darrieus 1 2 

for Giromill 3 4 

Figure 8 : Angle of attack vs. position 

for Darrieus 1 7 

for Giromill 3 4 

Figure 9 : Shaft alternatives 

9a : Stationary Inner Shaft 3 7 

9b : Stationary Outer Shaft or Guy wires 3 8 

9c : Single Rotating Shaft 3 9 

Figure 10 : Tower Shadow Effect 4 0 

Figure 1 1 : Darrieus Machine Configuration 5 0 

Figure 12 : Giromill Machine Configuration 5 1 


Figure 13 : Deployment Positions 

13a : Storage Positions 5 2 

13b : Blade Connectors 5 3 

13c : Details and Dimensions 5 4 

13d : Giromill Connections 5 5 

Figure 14 : Accurate Troposkien Shape 5 6 

Figure 15 : Pressure Distribution on Blades T54 

Figure 16 : Darrieus Airfoil T62 

Figure 17 : System Electronics 5 7 



k 

k 

EL 

L 

L 

l 


Statement of Purpost 


The purpose of the Mars Vertical Axis Wind Turbine group was 
to design a vertical axis wind turbine which would provide power to 
an unmanned Mars research station. The research package for such 
a station might include meteorology equipment, a seismometer, 
surface chemistry equipment, and imaging equipment. The power 
requirement of this package is one watt while in normal operation, 
non-transmitting mode (p. A-l to A-2). It is assumed that a lander 
is provided as a platform, with landing characteristics simi ar to 
those of the Martian Egg Lander EM 569 project of 1990 [19], and 
that the lander is equipped with a battery-type energy storage 
device The turbine should be able to operate in the Martian 
environment for three years, which is the amount of time between 
launch windows to Mars; continuous data is desirable, and it would 
be three years before another research package could be landed on 
the surface. A complete list of the needs and functions of the design 
as well as a FAST diagram appear on pg. A-3 to A-5. 


Environmental Concerns 


The Martian atmosphere differs greatly from the Earth's 
atmosphere. These differences strongly affect many aspects of our 
design. The most important differences are the temperature, 
pressure, and density of the air and the surface gravity. These 
values for Mars are: 


Temperature: 
Pressure: 
Density of air: 
Surface gravity: 


140 - 240° K 
6.5 - 10.5 millibars 
1.667 x 10' 5 g/cm 3 
3.70 m/s^ 


The values for pressure and density on Mars are between 1/75 
and 1/100 of their respective values on Earth. Even at these low 
densities, it has been shown that the Martian atmosphere may be 
modelled as a continuum rather than as a collection of discrete 
particles [5]. This fact enables us to use the same laws of 


aerodynamics as we would use for a wind turbine on Earth. J -j J— 

Another important environmental concern is that of dust Jk 
storms on Mars. The low density Martian air is able to transport r a 

only small dust particles (approximately 0.1 mm). The dust storms \ 
affect our design in two major ways. First, the airfoils will have to be 
made of relatively hard material to resist abrasion. -^Second, a shield 
will be used to protect the generator and the gears. 

The final environmental consideration is that of win ibeeds. 


The next section covers this topic in detail. 
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3¥ind Speed Design Decisions 

Two important decisions needed to be made regarding the 
Martian wind speeds. First, a range of wind speeds likely to occur 
most of the time needed to be determined. A good estimate for the 
maximum wind speed likely to occur was also necessary. 

The range of wind speeds needed to be determined in order to 
optimize the aerodynamic performance of our wind turbine. Because 
there is little data on Martian wind speeds, this decision needed to be 
based on a combination of analysis of the data and engineering 
judgement. The existing data consists of measurements taken at the 
Viking lander site and several meteorological studies. 

It is quite difficult to make any decisions based on the Viking 
data. The experiments measured north-south winds separately from 
east-west winds with no correlation between the two. The best 
estimate that we could obtain from this data is that the wind speed 
was greater than 4 m/s about 80% of the time [see pg. A-6], 

Studies done by meteorologists [5] show that the Viking lander 
site was a non-optimal site as far as wind speeds go. They are able 
to locate several areas where they believe the average wind speed to 
be above 8 m/s, and a couple where the average could be as high as 
14 m/s [5]. They also state that the global average wind speed is 6.5 
m/s. 

Since the mission will have other purposes other than just 
generating power from the wind, it is doubtful that the landing site 
would be at the optimal wind speed location. However, we believe if 
wind energy is to be used, the landing site should be a place of at 
least moderate wind speeds. Therefore we decided to design for an 
intermediate value between the Viking data and the meteorological 
studies. 

The actual numbers that we have selected are listed below 

1) Generate 1 Watt of power from all wind speeds 

greater than 6 m/s ^ ' 

2) Operate at peak efficiency for wind speeds between 

5 m/s and 8 m/s / 




3) Generate some power for all wind speeds greater 
than 4 m/s y 

We believe that these numbers will provide the best overall 
performance of our wind turbine. We see that even if we landed at a 
site similar to the Viking site, we would still be generating some 
power 80% of the time. 

The maximum wind speed likely to occur is very important for 
determining the extent of structural support required by our turbine. 
Wind speeds of up to 100 km/hr have been recorded on Mars for 
only a few hours [ 3 ]• Meteorologic studies estimate the maximum 
wind speed likely to occur on Mars at 250 km/hr [ 4]. 

Once again, engineering judgement needed to be combined with the 

data. 

The meteorological study here is somewhat questionable due to 
some of the assumptions it has made [ 41 - Also, 
expected to operate for only three years, so the 
seeing wind speeds of 250 km/hr is quite small. We decided to 
design our wind turbine to withstand all wind speeds up to 150 
km/hr. 




Selection of Wind Turbines 


Vertical-axis wind turbines appear to be the most feasible 
means for generating power from Martian winds. Since vertical-axis 
wind turbines are always oriented into the wind, there is no need for 
vanes to rotate the blades into the wind. This eliminates the 
presence of large gyroscopic forces found in horizontal-axis wind 
turbines. 

Several vertical-axis wind turbines were considered. The five 
possibilities that were given the most thought were: 

1) O-Darrieus 

2) A-Darrieus 

3) Giromill 

4) Split-Savonius 




A turbine that is powered by aerodynamic drag forces, such as 
a Split-Savonius, would need to be very large to extract enough 
power from the thin Martian air. Therefore, this choice was 
determined to be less desirable than the other four. 



Of the three lift-type turbines, the O-Darrieus and the Giromill 
were determined to be the best choices. Also, in order to minimize 
mass, each turbine will have two blades. Both the <I>-Darrieus and 
the Giromill possess strong advantages that were not found in the A- 
Darrieus. The main advantages of the two turbines are listed below. 


<P-Darrieus 

1. The blades are designed 

to eliminate bending stress 

2. Light weight blades can be 
built since no bending 

stresses are present 


Giro mi ll 

1. Blade material positioned 
for maximum aerodynamic 
torque 

2. Has deployment benefits 
since relatively easy to 

"compress" into small area 


A combination of the Split-Savonius with one of these two 
machines was also considered. This type of design would enable the 
turbine to be self-starting. However the large size of Split-Savonius 
that would be required prevented us from selecting this option. The 
use of the generator to start our turbine will enable a much lighter / 
design than if a Split-Savonius were used for this purpose. S & ^ ( 
For convenience, the rest of the report will drop the <|> prefix 
from the <}>-Darrieus. 
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Darrieus Blade Design 

• Size Determination 

• Development of troposkien shape for 
varying cross-sectional area of blade 

• Optimization of Shape 

• Airfoil Selection 

• Material Selection 


Structural Analysis 


Introduction 


When G. J. M. Darrieus received his patent for a vertical-axis 
wind turbine [2], he stated that the blades should take "the form of 
a skipping rope." A blade placed in this shape will be in a state of 
nearly pure tension, i.e. negligible bending stress. By eliminating 
the bending stress, the amount of material required for support is 
drastically reduced, and very light-weight blades may be designed. 
The word troposkien (from the Greek: xpojios, turning and 
oxoiviov, rope) is now used to describe the shape assumed by 
perfectly flexible cable that is spun about a vertical axis at constant 
angular velocity. 

All previous Darrieus turbines have been designed with a 
constant cross-sectional area along the length of the blade. 

Although these solutions ease manufacturing of the blades, much 
material is wasted. From previous designs, it has been found that . \ ? 
the tension, and therefore the stress, varies with horizontal ^ 

position squared [1]; the stress is highest at the top and lowest in 
the center. It was noticed that material could be saved by varying 
the cross-section along the length of the blade. Namely, it would 
be desirable to place more material where forces are higher, and 
remove some material where forces are lower. Our goal was to 
have a constant stress everywhere along the length of the blade, 
thus eliminating wasted material. Since the shape and the stresses 
depend upon how we vary our cross-sectional area, we realized 
that we could probably could not achieve exactly constant stress, 
but would try to come as close as possible. 





Size and Tin-Speed Calculations 


Using wind energy theory, together with experimental results 
from existing Darrieus turbines (see pg. T-l to T-8), we were able 
to determine the necessary size for our Darrieus turbine. In order 
to generate 1 Watt from a 6 m/s wind, the blades must have a 
swept area of 2 m 2 . Due to the low density of the Martian air, this 
is about 100 times larger than would be required on earth. The 
following figure shows what is meant by area swept. 

Area Swept , 


The aerodynamic performance of a wind turbine is strongly 
effected by its tip-speed. The tip-speed of a windmill is the speed 
at which the blades are moving ^ Since^ for a Darrieus blade, 
every position has a different speed, the tip-speed is defined as the 
maximum speed of the blade. This is by definition, the product of 
the distance times the angular velocity a (the distance b can be 
seen on Figure 0, pg. T-10). Closely related is the tip-speed ratio, 
which is the tip-speed divided by the ambient wind velocity. 

Using a suggested wind-speed probability relationship, the 
tip-speed was determined so that our turbine would generate the 
maximum average power (see pg. T-l to T-8). The average power 
was maximized so that we are able to recharge our batteries with 
the most energy over tirneQ^rhis was done while also assuring the 
generation of 1 Watt from a 6 m/s wind speed. The required tip- 
speed was found to be bco=31.3 m/s. 
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Troposkien Shape and Optimization 


The equations governing the shape of a varying cross- 
sectional troposkien are developed in complete detail in the 
appendix (see pg. T-9 to T-17). After considerable manipulation, 
the equations may be reduced to one first-order differential 
equation, one unknown, and two boundary conditions. The 
equations are in the form shown below. 


^ = f(r.b,C,K) 


r(z=a) = 0 


dr 

dz 


(z=0) = 0 


where 


Z 



r = horizontal coordinate 
z = vertical coordinate 
a = height, see above figure 
b = width, see above figure 

£ = a measure of how much we vary our cross-section along 
the length of the blade 

• a low C (C“0 means that there is a strong difference 

in cross-sectional area between the top and the 
middle of the blades 

• a high £ (£=20) means that the cross-sectional area is 

nearly constant along the length of the blade (this 
value for £ will be used to approximate the constant 
cross-sectional area solution) 

K = a parameter which depends on many terms, the most 
physically significant are: 


K is propof-KoAoJ 4-o aflju Ia*" veiocty 



• K is proportional to mass per unit length of blade 
f — a function, the actual function may be found on pg. T-17 


The values for a and b are not independent, they must be 
chosen so that the area swept is 2 m 2 . In order to solve for the 
shape r(z), we are free to choose any combination of a and b that 
sweeps out 2 m 2 , and any value of £. Since we have this freedom, 
we must select the values that give the best results. 

What our solution needs to minimize is the required 
arclength of the blades so that the area swept is 2 m 2 . The 
arclength is a measure of how much material is needed and 
therefore minimizing the arclength will minimize the blade mass. / 

The function f in the differential equation is quite complex 
(see pg. T-17) and no analytical solution exists. Therefore, we have 
developed a numerical solution with a computer program.*^ Our 
computer solution will be developed for several values of a and b. 
The selected values of a and b will not sweep out an area of exactly 
2 m 2 . Therefore our attempt will be to maximize the value of area 
swept divided by arclength; this is essentially a power-to-mass 
ratio. This concept can^learly be seen by observing the tabular 
data (pg. T-28 to T-31). Once this value is optimized, we may 
adjust a and b slightly in order to obtain the swept area equal to 2 
m 2 . 

The computer solution (pg. T-18 to T-27) is carried out by 
first selecting values for a,b, and £, and then solving for the shape 
r(z) and the parameter K. The shape is determined pointwise and 
may be developed for as many points as desired. Since the 
equation is a first order differential equation with two boundary 
conditions, we need to leave K as an unknown in order to 
guarantee that the solution will "fit" both boundary conditions. 


The program uses a combination of the fourth order Runge- 
Kutta method for numerical integration and the bisection method 
for root finding (see pg. T-24 to T-27). The program is able to 
determine the shape, r(z), and the value of K simultaneously. For 
each selected value of the solution is repeated for 30 different 
combinations of a and b that sweep out nearly 2 m 2 . The whole 
procedure is then repeated for different values of with the goal 
being a maximum ratio of area swept to arclength. * 

After trying several values of £, we notice that the lower the 
C, the higher the value of area swept to arclength for all values of a 
and b. The following graph shows this quite clearly. The ratio of 
area swept to arclength is shown on the vertical axis. The height- 
to-width ratio (a/b) is plotted on the horizontal axis. The curve for 
C=20 is essentially the case of a constant cross-sectional area along 
the length of the blade. 
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From these curves we see that by varying the cross-section, 
we get the added benefit of a better shape. This makes our results 
better than we had even hoped for. The concept of obtaining a 
better shape also makes sense physically. By placing more 
material near the top, we increase the centrifugal force there, and 
obtain a shape that is more square-like. This shows up well on the 
following graph. The graph shows the shape obtained by varying 
the cross-section (£=1.3), the shape obtained with a constant cross- 
section (£=20), and a circle for comparison. As expected, the 
varying cross-section case is "pushed out" more near the top of the 
blade, and the constant cross-section case is "pushed out" more 
towards the center of the blade (which is the bottom on this 
graph). All of these curves have the same arclength. 

COMPARISON 


Scale: Legend: 
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We see from figure 2 that for each value of £, there is an 
optimal combination of a and b. We are then left with determining 
the best value for £. By definition, £ has a minimum possible value 
of approximately 1.0. With this fact and the above curves in mind, 
we see that we want £ very close to 1.0. 

Besides having the best shape, we must remember our 
original intent to have nearly constant strength (i.e. constant 
stress) along the length of the blade. It is important to understand 
the differences between the two measures. A constant strength 
reduces required material along the length of the blade, whereas a 
better shape reduces the length of the blade. 

After several iterations, we see that if f; is too close to 1.0, we 
are not able to approach constant strength and that if £ is much 
greater than 1.0 we quickly lose the better shape and the ability to 
obtain constant strength. The optimal value, determined by trial 
and error, was found to be £= 1.3. 

With this value selected, and the optimal values for a and b 
determined from the tabular data (pg. T-28 to T-31), we have the 
complete solution for our shape and the value of K. With b known, 
we are able to determine the required angular velocity from our 
earlier specified tip-speed of ba>=31.3 m/s. Then with K and oo 
known, relationships may be written in order to determine the 
tension and stress at any point along the blade. This is all done on 
pages T-32 to T-36. 

From these relationships, we are able to plot the tension and 
stress along the length of the blade. The values are normalized for 
ease of understanding. The tensile value at any point is divided by 
the tension at the blade's midpoint (z=0), which is the minimum 
tension. The stress value is also divided by the stress at the 
blade's midpoint. These ratios are written as T/T 0 and o/a 0 , 
respectively, where the subscript o refers to the value at z=0. 

These ratios are measured along the vertical axis on the following 
plot; Along the horizontal axis is the non-dimensionalized vertical 
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coordinate (z/a). We see that we have taken a maximum loading 
ratio (T/T 0 ) of nearly 4 and reduced it to a stress ratio (a/a 0 ) of 
approximately 1.5. We also see that we have come very close to 
achieving constant strength (o/a 0 =l) throughout the length of the 
blade. For the traditional, constant cross-sectional area Darrieus 
designs, the stress ratio would follow the same curve as the tension 
ratio. 


TENSION AND STRESS RATIOS 

^ ca ^* : Ltftni; 

H. lack - • - TENSION' RATIO 

V. lack • 1.706* STRESS RATIO 
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Although all of these figures give great promise to the idea of 
varying the cross-sectional area along the length of the blade, we 
needed to estimate the total saving of mass to see if the idea’ was 
really worthwhile. Remember that benefits occur for three 
different reasons. First, we reduce pia^s by varying the cross- 
section along the length of the blade '"and removing material where 
it is not necessary. Second, we obtain a better shape (shorter 
length of blade) which requires less material.^ Finally, by reducing 
the total blade mass, we reduce stresses everywhere.^ When our 
solution was compared to the constant cross-sectional area solution 
we found that by varying the cross-sectional area we reduced 
blade mass by over 34% while also achieving a 54% reduction in Q 
maximum stress (see pg. T-36). These results verify that it is 
worthwhile to vary the cross-section along the length of the blade. 



Darrieus Airfoil Selection 
Chord Leneth 


The optimum performance of the wind turbine depends upon the size and 
shape of the machine itself and the size and type of airfoil used for the 
blades. A measure of the turbine’s performance can be given by 


Cp = 0.25n(c/r)kXV2 . 0.5n(c/r)CdX^ 

where: 

Cp = coefficient of performance, a measure of efficiency 
n = number of blades for the turbine 
c = chord length of the blades 

r = distance from the chord line to the center line of rotation 
X ■ the tip-speed of the machine 
V = the incident wind velocity acting on the machine 
Cd = the average drag coefficient for the blades 

The incident velocity, V (Fig. 6 and appendix p. T-41; [24]), is given by the 
equation 


V = 1 - 0.0625n(c/r)X(k+3Cd) 


Substituting this equation into that of Cp, differentiating with respect to 
the tip speed, X, and assuming the atmosphere acts as an ideal fluid (Cd=0) 
results in the expressions 


and 


X - 16r/(3nck) 
V = 2/3 


The latter simply means that during optimal operation, the incident 
velocity acting on the wind turbine is equal to 2/3 the free stream 
velocity. For the darrieus machine, the first equation, used to find the 
chord length, is now a function of two variables, c and r; unlike the 
giromill. Considering that the torque generated by the the outermost part 
of the blade will be much greater, than that of inner sections, the value of r 


- 'S'- 


r 



used was 0.839 m. Substituting values for X.r.n, and k results in an 
optimum chord length of 6.83 cm. 



o 10 20 30 40 


■ cp 
• cp 

■ cp 
e cp 

■ cp 


Free stream velocity 


"TKd, p re,<i\'c4-e A Curve cocoas eXpc* ;as4*\4^1 

resu.1^ Aone, on | 'Derr.'ew.^ +wrkin£.S, An Id 

<>f -4es>e e^j>eri*w4.n+A.l curves «>*v - ~T~4. 


To 


,art 









Thickness 


The thickness of the airfoil can be determined by finding the maximum 
value of the average coefficient of moments (Fig. 7) which is given by 

Cma = (c/r)[(p/2)V2-CdX2] 


By inspection, it can be seen that the thickness of the airfoil should only be 
dependent on finding the lowest value of Cd- As before, the radius also 
affects these calculations. However, as r decreases, the contribution to the 
moment decreases as well and so we make the same assumption of r as 
when finding the chord length. Using experimental data [8] to compute a 
value for these drag coefficients gives the following as examples: 

Airfoil type Qj 


NACA0006 

NACA0009 

NACA0012 

NACA1410 


0.0098 
0.0079 
0.0076 ^ 
0.0081 


As can be seen, there is very little difference among these drag values 
which results in moment coefficients which differ by only thousandths; the 
NACA 0012 airfoil being only slightly better than the others. Initially it 
was thought ignoring these differences and using the thinnest section 
possible (to cut mass requirements) would be best. The range of usable 
angle of attacks (Fig. 8 and appendix p. T-42; [24]) for each airfoil, 
however, does show that the NACA 0012 section is the best choice, 
especially at higher wind speeds. 


Lift and Drag 

Due to the extremely low density of Mars' atmosphere and the small 
planform area of the blades (0.18 m^) the lift and drag forces acting on the 
blades are small. At a wind speed of 6 m/s the lift force on the blades is 
only 0.72 N (appendix p. T-42). The drag forces have been neglected due 
to the fact that they are typically on the order of one percent of the total 
lift force. 
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Darrieus Material Selection 






/ ' LS' [ 

The material selection for the blades was limited % twb, mam 
factors, one being the low temperature range (140K-240K), and the 
second being the need for a high strength-to-weight ratio. The 
materials collected within those limits are found in Table 1 and Table 

2 on pg. M-l and M-2 in the appendix. 

The Darrieus blades are in pure tension while the machine is 
operating, thus a light weight material with strong tensile properties 
is needed. These two criteria are best exhibited by unidirectional 


1 o iivvuvv* • * ' 

composites ( Table 2 (pg. M-2) in the appendix shows properties of 




° 
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such composites). ^ - C T ^ 

Kevlar had the highest strength-to-weight ratio of all the 
materials. It also exhibits superior properties of fatigue. KevlarJ^ 1 
suited the need for the Darrieus blades and thus was chosen for the 

bfedTHaterial. Verification of the material selection can be found on 

{S r* ( S 

oaee M-4 in the appendix. / > > y a ( d u 

All connector pins are made out of titanium (alloy 6A14V. This 
material was chosen because of its extremely high shear strength. Its 
low coefficient of expansion is also very important so that it won’t 
interfere with the materials it is pinning. Properties of this titanium 
alloy can be found in Table 1 (pg. M-l) in the appendix. 
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The solution of the troposkien equation for a given set of 
parameters was completed. These parameters constitute a set of 
"governing equations" that determine the dimensions of the blade. 
The "governing equations" have been included in the appendix(T-33, 
T-34). The subsequent analysis of the blades included the following: 

(1) Centrifugal loading analysis, 

(2) Aerodynamic loading analysis, 

(3) Thermodynamic loading analysis. 

These will be discussed further in the following, additionally, sample 
calculations have been included(T-43 - T-62). The troposkien blade 
specifications will be discussed first. 

Troposkien blade specifications 


The final configuration of the troposkien blade was a hollow 
shell with approximately constant wall thickness. Many alternative 
configurations were considered(T-60,T-61), but in the end the thin, 
hollow shell proved the best choice(T-63). 


Outer shell shape 


- NACA0009 







At the equator 

cord = 7.4 cm 

wall thickness = 0.5 mm 


At r = 6 cm 

cord = 9.0 cm 

wall thickness - 1.0 mm 


total mass of each blade = 0.486 kg. 

Addition information on the geometry of the blade can be found in 
the appendix(T-43 - T-45). 
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Centrifugal Loads 


The dominant loads acting on the troposkien airfoil are the 
centrifugal forces. The troposkien blades are in a state of pure 
tension for centrifugal loading. The solution of the troposkien 
solution yields specific parameters which must be satisfied; these 
parameters may be found on pg. T-33 and T-34. Manipulation of 
these equations yields specific information about the stress in the 
blade at any location (T-39). Remarkably, we find that the tensile 
stress is completely independent of the cross-sectional area, as long 
as the area is varied according to the mass per unit length equations 
on pg. T-14. This is because as the cross-sectional area is increased, 
the centrifugal force is increased, and the tensile force increases. For 
example, the uniaxial stress acting at the blade's equator is 


G 0 = 61.1 * 8 — T 0 / (A cs ) 0 

where: 

T 0 = tension at equator 

8 = mass density (mass/length^) of material 

(Acs)o = cross-sectional area at the equator 


The density of the blade material is 1380 kg/m^. This yields a 
uniaxial stress of 84.3 kPa at the equator. The maximum stress is 
129.8 kPa, occurring at the roots of the blades. Because of the high 
strength of Kevlar, the factor of safety under normal operating 
conditions is 1 0,870 ^ The only way to reduce the factor of safety is 

to use a weaker and/or heavier material! 

The source of greatest concern with respect (rijp the centrifugal 
loads is the location of the pin that will secure the blade to the hub of 
the shaft (T-46 to T-48). A stress concentration will occur around j (j 
the hole. The stress concentration factor [20] is defined as, L# 

k = 2 - (1-w/d) + [1.5(l-w/d) / (1+w/d)] c/ 

The stress intensity factor evaluated at the pin location is ^k=l,725 j/ 

The maximum stress is computed by G max = k P/t(w-d). The 
load P is determined to be T/2 = 9.0 N, so we obtain 
O max = 181.1 MPa 


' z 


The ultimate stress for the material is given to be 1400 MPa, giving a factor 
of safety of 7.7. 


Aerodynamic loading 

The aerodynamic loads acting on the airfoil were found to be 
negligible compared to the centrifugal loads. Considered in the analysis, 
were the pressure acting over the surface of the airfoil, which could 
potentially deform the airfoil shape, and the bending moment acting on the 
blade due to the total aerodynamic load. 

The pressure acting on the surface of the blade was analyzed by 
modeling the blade surfaces as infinite plates fixed on the leading and 
trailing edges. 


The equation for the maximum deflection of the center of the plate[22]is 

ymax = a q b 4 / E t 3 

where 

a = 0.0285 

q = pressure acting on the surface 
b = cord of airfoil ^ 0 ? j '- j 

E = flexural ^£ odui es)of material 
t = wall thickness of the blade 

Assuming the maximum deflection of the plate to be equivalent to 
half the maximum thickness of the airfoil, a pressure can be calculated. 
This pressure acts as a residual force against any restoring pressure. 
Thus, for ymax = 0.00333 m, 

q = 38.96 MPa 
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The pressure acting on the blade due to the aerodyn ami ra of the 
airfoil is estimated(T-51 - 1-55)[8] and subtracted from q. 


Q = q - p = 38 94 MPa 


This yields a maximum deflection of 

ymax2 = 0.003328 m 

The difference between the maximum deflections is 0.05%. 

The aerodynamic loading was estimated to be 0.06 N-m about the 
vertical axis. The moment acting on the blade at r = 0.06m would be 
approximately equal to that acting about the center(T-56 - T-57). 

The relevant moment of inertia and centroid at r = 0.06m is 

Iyy = 514.08(10)-9 m 4 

x = 0.042c 

Substituting into a = Mx / Iyy, yields 

a = 6.1 kPa. 

This stress is less than 5% of the uniform stress due to the 
centrifugal loads. 


Thermody namic loading 

The martian atmosphere experiences a diurnal temperature 
change of -100 K. The strains that result from these temperature 
fluctuations is negligible(T-58, T-59). 

The material has an axial and a transverse coefficient of 
expansion. They are: 


E axial = 100 K (-2.0(10)-6 / K ) 
E trans = 100 K (60(10) - 6 / K ) 


In the axial direction, the total length is 2.54 m. In the transverse 
direction the length is approximately the cord length. The total 
change in overall length is 


6-53.4 ura, axially; and 
6 = 444 (im, transversely. 




r if’ 




I 




The expansion in the axial direction is negligible, and the expansion 
in the transverse direction is less than 1%. 
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Giromill Blade Design 


• Size Determination 

• Strut Positioning 

• Blade Structure 

• Strut Connections 

• Stresses 

• Airfoil Selection 


Material Selection 


CIROMILL 


Dimension ? and Tip Speed Rati g 


In determining Giromill dimensions the initial step comprised of 
an approximate Swept Area(Asw) (see Pg.Gl) needed to produce one 
Watt of mechanical power. From Asw dimensions of height an 
radius of Giromill were obtained. 

From wind energy theory, it is possible to deduce a rough idea of 
Ajw. By considering kinetic energy of moving air, the following 
equation for maximum mechanical power(Pioech) that can be 
extracted was derived:(for complete derivation see pg.T-S) 

Eq.(l) Pmech=Cp(.5*p*Asw*Vo A 3) 


Cp=>coefficient of power 
p=>density of Martian air 

Vo=>free stream velocity on Martian surface 


Wanting to maximize Pmech, a peak Cp and related Tip Speed 
Ratio(T) where needed. Numerous numerical data showed an 
approximate Cp/m.i=0.5 in accordance with a T-3.0. [9](see pg. uz; 

Solving Eq.(l) for Asw, setting Pmech=1.0 Watt, and plugging in 
appropriate quantities gives: 

Eq.(2) Asw=1.389 [m A 2] (see pg.G3 for calculation) 

It is also know that the Asw can also be represented as a function 
of height(h) and radius(R) of the Giromill Eq^3) below Another 
equation which is a function of h and R is the perimeter of the 
Giromill(P) Eq.(4) below. 


Eq.(3) Asw(h,R)-2hR-1.389[m A 2] 

and 

Eq.(4) P(h,R)=4R+2h 


If Eq.(4) is minimized to ensure smallest amount of material 
used, Eqs. (3) and (4) can be solved for Giromill dimensions h and R. 


R=.5546[m] 

h=1.1786[m] 
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(see pgs.G3-G4 for calculation) 

The T related to Cp/max can be used to obtain the operating 
angular velocity of the Giromill, with the correct R determined 
above. The equation for T is as follows: 

Eq.(5) T= (i)R/Vo 

o)=> operating angular velocity of the Giromill 
Setting Eq..(5) equal to 3.0 and solving for to: 

co=32.45[r/s] 

(see pg.G5 for calculation) 

These first calculations for dimensions and angular velocity 
became very important for later calculations in Giromill design. 


Strut Positioning 






The decision for a two strut per blade desijgn^w^^determined 
from two factors. Factor one was two struts 'dre" a simpler design 
/thieipthree or more. Factor two was two struts reduced the 
maximum bending moment in the one strut case by 17. 1%. (see 
appendix pg.G15) 


For the two strut design the positions of the struts along the 
blades was crucial. The position was in-direct relation -to-thfr bending 
moments and stress on the blade. •""’ J ^ 


In the blade and strut arrangement, the blade was modelled as a 
beam with a uniform distributed load(w) acting in the plane of the 
struts. (see appendix pg.G9) Looking at the bending moment diagram 
we see three large moments: 


MA=MB=wa A 2/2 

MC=(wh/8)(-h+4a) 

a=>position of struts from end of blade 
h=>length of blade 

(see appendix pg.G12-G13) 
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With the maximum bending moments determined, it was 
desirable to minimize them as much as possible. This in turn 
minimizes the stress on the blade. If the moments MA and Me can 
be balanced the desired minimal stress will result. By determining 
the ratio of a/h to make the following equality true: 

Eq.(l) |MA|=|MC| 

the smallest stresses possible are achieved. 

Plugging the appropriate equation for the moments in Eq.(l) from 
above The following equation for a and h result: 

Eq.(2) a A 2+ha-.25h A 2=0 

Solving Eq.(2) give the ratio a/h=. 207 or a is 20.7% the length of 
the blade under a uniform load. (see pg.G13-G14) 


Blade Structure 

A NACA-0012 airfoil shape was chosen for reasons to^be 
mentioned later under the title Giromill Airfoil Selection. The 
material aluminum boron was chosen for the Giromill, which includes 
the blades, for reasons which will be mentioned under the title of 
Giromill Material Selection. Z' 

When considering blade structure, centrifugal loads acting under 
operating conditions become the key design consideration. The 
aerodynamic loads are negligible compared to the centrifugal 
loads. (see pg.G18) 

Centrifugal distributed load: 

Eq.(l) W=(CO A 2R/h)Mass|one blade 
(see appendix pg.G19) 

This relates to the maximum bending moment on the blade of: 

Eq.(2) MA=wa A 2/2=((<o A 2Ra A 2)/h)Mass|one blade 
(see pg.G12) y 
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Notice Eq.(2) Ma is a function of mass, if this mass can be reduced, 
the stress on blade will also reduce. 

To reduce the mass the blade was hollowed out to a thickness of 
.5[mm]. This dimension was determined on a manufacturability of 
aluminum boron. 

The consideration of a ribbed blade was omitted due to the fact 
that^ with the rib^ stress and deformation of blade increased. 

Stresses fusing ax on pg.G21) 
w/o rib (Jx|w/o=1.6297* 10 A 7[Pa] 

w/rib Ox|w/= 1 .6899* 10 A 7[Pa] 

Deflectionf see calculation on pgs.G28-G29) 

yE|w/o=(.9644)yE|w/ 

The final overall dimension of the blade are as follows:(see pg.G4, 
pg.G30 and section Giromill Airfoil Selection) 

C(cord length)=8.13[cm] 

T(Max. thickness)=9.756[mm] 

D(blade thickness)=.5[mm] 
h (blade length)=l. 17857[m] 


Connections of Struts 

To connect the blades to the struts three points must be consider. 
Point one the struts must be pinned for deployment and to reduce 
bending moments on blade/Point two struts must be positioned at 
42.04% of the cord length from the leading edge. This will insure no 
resulting moment due to > centrifugal loading. Point three struts must 
be positioned at 20.7% of blade length from end of blade.(see pg.G14) 

To comply with point one small cylinders of radius=1.5[mm] and 
length=7.0[mm] are attached to the end of each strut. This small 
cylinder is then fitted into a U-shaped fitting, which is welded to the 
blade. See appendix pgs.G35-G37 for complete dimensions and shape 
of pin connections. 
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Position for points two and three are found in appendix pg.G34. 


Stresses and Factor of Safety 

In the final stress analysis there are two points to consider. Point 
one is the the centrifugal loads dominate the aerodynamic loads. (see 
pg.G18). Point two the centrifugal loads act perpendicular to the 
blade and at the center of mass. 

Keeping these two points in mind, the maximum stress(ax) found 
in the blade will occur at the point of maximum bending moment 
Calculated ax during operating condition: 

Ox=l. 62960792* 10 A 7[Pa]= 16. 296[MPa] 
with a factor of safety of: 


F.S.=81.0 


(see pg.G41 for sample calculation) 




Checking shear stress in the blade is approximately during 
( operating conditions: 

x=4.1[MPa] 

with a factor of safety of: 


F.S =19.5 


(see pgs.G42-G44 for sample calculation) 

Notice the large factor of safety. Due to the nature of the 
centrifugal loads and moment of inertia dependance on blade 
thickness(D) this can not be helped. The stresses depend on (D) in 
such a way that they decrease with decreasing (D). (see pg.G39 for 
justification) 

The final consideration looked at concerning the blade was the 
pressure load due to aerodynamic loads. It was thought that with 
such a small thickness of the blade (D=.5[mm]), this pressure load 
may deform the blade. This was modelled and determined to be 
negligible. (see pg.T51-T55 for calculations) 
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The stresses concerning the struts and connections are of very 
high factors of safety. This is due to the small forces applied to the 
strong aluminum boron material. See appendix pg.G31-G34 for strut 
considerations and dimensions. 

Total Mass of Blade s and Striits 

The Giromill contains two blades, four struts, and four pin 
connections. The resulting mass of Giromill(MT): 

MT=2MB+4Ms+8MU+4MCL=0.42887[Kg] 

(see pgs.G45-G46) 


nirnmill A irfnil Selection 
Chord Length 


The optimum performance of the wind turbine depends upon 1 e or^he^ 

shape of the machine itself and the size and type of airfoil used for 
blades. A measure of the turbine’s performance can be given by 


C p = 0.25n(c/r)k\v2 . 0.5n(c/r)CdX 3 


where: 

C p = coefficient of performance, a measure of efficiency 
n = number of blades for the turbine 

c = chord length of the blades . 

r = distance from the chord line to the center line of rotation 

X = the tip-speed of the machine 

V = the incident wind velocity acting on the machine 
Cd = the average drag coefficient for the blades 

The incident velocity, V (Fig. 6 and appendix p. G-54; [24]), is given by 
the equation 


V = 1 - 0.0625n(c/r)X(k+3Cd) 


Substituting this equation into that of Cp, differentiating with respect to 
the tip speed, X, and assuming the atmosphere acts as an ideal fluid (Cd ) 

results in the expressions 


and 


■ 16r/(3nck) 
V = 2/3 


/ 


The latter simply means that during optimal operation, the incident 
velocity acting on the wind turbine is equal to 2/3 the free stream 
velocity. For optimum performance, the first of these results as on y 
variable, c. For the case of the giromill. Substituting values forX.r.n, and 
results in an optimum chord length of 8.13 cm. 


- 31 - 



Thickness 


The thickness of the airfoil can be determined by finding the maximum 
value of the average coefficient of moments (Fig. 7 and appendix p. G-54 to 
G-55; [ibid.]) which is given by 

Cma = (c/r)[(p/2)V2-CdX 2 l 

By inspection, it can be seen that the thickness of the airfoil should only be 
dependent on finding the lowest value of C d - Using experimental data [8] 
to compute a value for these drag coefficients gives the following as 
examples: 


Airfoil tvoe 

Cd 

NACA0006 

0.0098 

NACA0009 

0.0079 

NACA0012 

0.0076 

NACA1410 

0.0081 

seen, there is 

very little difference 


/) 

5 ^ 
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which results in moment coerncienis which umci 

NACA 0012 airfoil being only slightly better than the others. Initially it 
was thought ignoring these differences and using the thinnest section 
possible (to cut mass requirements) would be best. The range of usab >e 
angle of attacks (Fig. 8 and appendix p. G-54 to G-55; [24]) for each airfoi 1, 
however, does show that the NACA 0012 section is the best choice, ^ 
especially at higher wind speeds. Peculiar to the darneus type wind 
machines is the fact that regardless of the wind speed, chord length, and 
thickness, there will always be a section of the blade that will be stalled 
during some portion of the rotation. As r decreases or as the free stream 
wind velocity increases, the tip speed X decreases. This results in an 
increasing angle of attack until the stalling angle of the airfoil is reached. 
At X = 1 the blade speed is equivalent to the wind speed and the lift and 
drag forces are equal to zero at that instant. At optimum levels, this point 
occurs at roughly 16 cm from the axis of rotation. 
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Due to the extremely low density of Mars' atmosphere and the small 
planform area of the blades (0.09m2) the lift and drag forces acting on the 
blades are small. At a wind speed of 6 m/s the lift force on the blades is 
only 0.23 N (appendix p. G-54 to G-55). The drag forces have been 
neglected due to the fact that they are typically on the order of one 
percent of the total lift force. 
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Giromill Material Selection 


The material selection for the Giromill was limited by two main 
factors, one being the low temperature range and the other being the 
need for a high strength-to-weight ratio. The materials collected 
within those limits can be found in Tables 1 and 2 on pg. M-l and 
M-2 in the appendix. 

The material selection for the Giromill blades and struts must also 
take in consideration bending, and shear. A high modulus of 
elasticity is also needed to withstand any deflection due to the 
location of the struts. The material which exhibited the best 
properties in all categories was boron reinforced aluminum. This is a 
metal matrix composite characterized by high pensile strength and 
shear modulus, dimensional stability, j oinability, high ductility, and 
toughness. For further verification of the Giromill material selection, 
see page M-6 in the appendix. 

All connector pins are made out of titanium alloy 6A14V. This 
material was chosen because of its extremely high shear strength. Its 
low coefficient of expansion is also very important so that it won't 
interfere with the materials it is pinning. Properties of this titanium 
alloy can be found in Table 1 on page M-l in the appendix. 


^ jvj£\<Ufu 1 1 rij J O 



Shaft Design 
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The loading on the shaft is as follows (pp. S-l to S-8): 

•Bending due to wind drag on the turbine blades 
•Torsion due to power transmission from blades to 
generator 

•Axial/buckling due to weight of blades ® 

The shaft configuration of Figure 9a was chosen because fatigue 
is eliminated from the shaft. The non-rotating inner shaft takes the 
bending moment, while the outer shaft takes the torsional load. 

Because the rotating shaft does not see any bending, a fatigue 
situation is avoided. 
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The bending moment of approximately 18 N-m led us to an 
aluminum-boron composite (A1B) inner shaft with a diameter of 

l 

10mm. The outer rotating shaft, also of AIB, has an inner radius of 
19mm and an outer radius of 19.5mm. The bearing races are 
integral to the shafts to avoid problems with adhesion over the low, 
wide temperature range. The dimensions of the bearings were based 
on an SKF ball bearing which was selected for its ability to handle 
small thrust loads [p. S-l 3]. The races are coated with a film of 
molybdenum disulfide for lubrication. 

A first approximation of the torsional natural frequency of the 
Dameus machine was found to be 2.5 Hz, while that of the giromill 
was 4.3 Hz (pp. S-9 to S-l 2). The main torsional forcing function for 
these machines is caused by "tower shadow," which is illustrated in 
Figure 10. Each blade passes through a region of reduced wind/ 
velocity, caused by the interference of the tower, once per 
revolution. When the blade passes through this tower shadow, it 
experiences a dip in lift and therefore in the torque transmitted to 
the shaft. Because there are two blades, the frequency of this dip is 
twice the rotational speed of the machine, or 5.6 Hz for the Darrieus 




Outer rotating shaft 



Inner stationary shaft 


Figure JL a: Stationary Inner Shaft Configuration 

Note that the inner shaft bends while the outer shaft • 
rotates and transmits power to the generator--no 
fatigue. 
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Figure _2_b: Stationary Outer Shaft or Guy Wires 

Note difficulty of power transmission from 
inner shaft to generator outside of outer shaft. 
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Tower Shadow 



Figure Id.:. Illustration of Tower Shadow. Our 
shaft is so small in diameter that this effect is 
minute at best. 
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and 5.2 Hz for the giromill. If the torsional natural frequency of the 
shaft is near the frequency of this forcing function, a damaging 
condition of resonance could occur; however, the amplitude of the 
forcing function depends directly on the width of the tower. Our 
tower is so slender (less than two percent of the total width of the 
machine) that it is questionable whether a turbulent shadow would 
even exist behind our shaft. Regardless, the resonant frequency of 
the Darrieus machine is considerably lower than the frequency of the 
forcing function; that of the giromill is closer but still less than half of 
the forcing function frequency. We have therefore concluded that 
the major vibrational forcing function of tower shadow should not 
present a problem. f 

The total mass of the shaft for the Darrieus is 0.396 kg, and for 
the giromill is 0.300 kg. Development of shaft loading and 
dimensions, and natural frequency calculations appear in greater 
detail in Appendix S. Figures 11 and 12 show the configurations of 
both machines. J . 
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Gears .Bearings and Lubrication 


/The bearings are integral to the rotating and stationary shafts. 

The power is transmitted from the rotating shaft to the generator 
through a gear system. The large gear is integral to the rotating 
shafu'lpid the small gear is connected to the generator. A gear ratio 
of \:6js needed to accommodate the generator. A detailed 
description of the gears and bearings is found on pages M7 and M8 
in the appendix. 

Lubrication becomes a complex problem when dealing with the 
low temperature range of Mars. Oils and greases do not operate at 
these low temperatures/ therefore solid lubricants must be used. Of 
the solid lubricants available. Molybdenum disulfide (M 0 S 2 ) is the 
best suited for our purposes. M 0 S 2 is effective in a vacuum, ■/ 
dependable at low temperatures, dimensionally stable, and is not 
damaged by radiation. Other solid lubricants , such as graphite, are 
ineffective in at least one of these area [l 1,1 2 L M 0 S 2 has a 
comparably low coefficient of friction and one of the highest wear 
lives (see Table 3, pg. M-3 in appendix). M 0 S 2 performs best when 
used as a film. To increase wear life, M 0 S 2 should be bonded to the 
surface with a phenolic resin. 

resin bonded M 0 S 2 film wear life= 9,860,000 (35ksi cycles) 

non-bonded M 0 S 2 film wear life= 103,680 

This M 0 S 2 -resin combination works well with aluminum, which is 
the material that we are using for our gears and bearings. The 
thickness of a typical film for our case would range from .07mm - 
,001mm. 


Generator and Electrical System 


The electrical system for this project has three primary functions 
These functions include: 

1. ) Rotational start-up ^ 

2. ) Speed control ^ 

3. ) Electrical power generation / 

To achieve each of these objectives simultaneously, a rather complex ^ 
system will be employed. A schematic of this system is shown in fig. 17.^ 
The system includes a miniature DC motor with a step-up gearhead, a 
speed control card with potentiometer, a battery bank with charging unit, 
an external anemometer, and two control devices. Each function, along 
with its individual components is examined in detail in the pages that 
follow. 


Rotational Start-Uft 

Neither the Giromill nor the Darrieus Vertical Axis Wind Machine 
(VAWM) are expected to be able to self-start with any degree of certainty. 
To overcome this problem, the motor will be used to start the rotation of 
both machines. 

The estimated wiiuK speed at which both machines would be able to 
produce energy is 4 m/s (see p. T-4). The anemometer would be used as a 
means to detect when the free stream wind speed reache<£ this cut-in 
value. Upon reaching the cut-in wind speed, the cut-in control would be 
used to forward-drive the motor and begin rotating the blades of the 
VAWM. 


To forward-drive the motor, the battery bank would have to be used 
as a power supply for the motor.-'iThis type of operation would present the 
potential problem of completely exhausting the batteries. To prevent this, 
another control device would be used that would limit the level to 
the batteries could be depleted,/^ 

' pi* ■ 

M r 



Speed Control 


Speed control is a very important part of our VAWM design. For the 
case of the Darrieus, the moment-free shape of the troposkien is only valid 
for one rotational speed. The maximum coefficient of power is also 
attained at a constant tip-speed-ratio, which, for a constant wind speed, 
requires a constant rotational speed. The final reason for constant 
rotational speed is structural integrity. If the machines were allowed to 
spin uncontrolled, the stresses developed would become self-destructing 
and eventually cause failure. 


To maintain a constant rotational speed, a commercially available 
speed controller card, the Instech 1100, will be used (see reference [25]). 
The rotational speed of the motor (back-driven to act as a generator) is 
linearly related to its output voltage by the generator's velocity constant. 
By measuring the output voltage and knowing the velocity constant of the 
generator, the card is able to use an external potentiometer to regulate the 
rotational speed of the generator to a specified value. 


Examining the figure on p. T-4, we see that output power increases 
with increasing wind speed up to a point after which the output power 
decreases. In general, we find that, "If the fixed speed load is able to 
accept the maximum possible mechanical power, no additional braking or 


loading is necessary as the wind speed increases above its^ rated vs^j 

By scaling the figure on p. T-4, we estimate a 0 

output capacity of 3.1 Watt^will be required for adequate braking. For S 




confirmation of these findings, ho\ 
performed. { 

Electrical Power Generation 


^ver, wind tunnel tests should be 
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For braking, we found that a 3.1 Watt generator is required. NASA 
also requires that the system operates at 12 Volts. This information allows 
us to examine available miniature motors. A commercial miniature motor 
catalogue indicates that a 12 Volt, 3.7 Watt motor would typically have an 
efficiency of 85%, a velocity constant of 714 r.p.m./volt, a no-load speed of 
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9,000 r.p.m., and a mass of 0.058 kg [26]. 

For the operating speed of 300 r.p.m. to generate 12 Volts, a gear 
ratio of 30:1 will have to be used. This will be achieved by a 6:1 increase 
from the blade rotational shaft to the shaft of a 5:1 gearhead. Again 
referring to a commercial motor catalogue, we fmd an efficiency of 80% 
and a mass of 0.065 kg for such a gearhead [26], ^ / 

Calculating the driving torque imposed on the central shaft by the 
blades, we fmd a value of 0.068 Nm (pp. P-1, P-2). At the normal ^ 
operating speed of 300 r.p.m., we find an output power of 1.48 Watts(p. P- 
2). This value does not take the efficiencies of the bearings into account, 
and will be slightly higher than the actual value. 


Other C onsiderations 


1. ) Conventional lubrication will probably not be effective at 

such extreme temperatures. Investigation of alternatives 
should be made. ^ 

2. ) Calculations show that at 15% inefficiency, a 3.7 Watt j 

motor would reach steady state 17.21 K warmer than / j 
the surroundings, (p. P-3) This does not present a / 
problem for overheating. / ' 

3. ) The total mass of the electrical system (excluding batteries 

and charging unit) is 0.160 kg. 
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Deployment 


The Darrieus and giromill wind machines will be stored in the 
lander in a very space efficient manner. 

DARRIEUS 

The Darrieus machine has blades that fold into the center shaft 
when stored. See figure 13a. for stored position. Upon landing on 
the surface of Mars, the lander lid will open and the machine will be 
elevated out of the lander compartment. The stationary shaft of the 
machine will be secured to the elevated platform which rises and 
acts as the base of the wind machine. This platform will have a blind 
hole that the windmill's stationary shaft will be inserted into. Once 
the elevator has deployed the machine from the lander 
compartment, the blades will fall into thei^ appropriate operational 
positions due to gravitational effects. Seemgure 11 for the elevated 
working position. The blades simply roll around a countersunk pin 
hole into their locked and working positions. This countersunk hole 
is higher at one side with an increasing slope leading to a slot that 
the connecting pin will lay in for operation. See 'figure 13b. for 
detailed drawing of blade connectors. The blades of the Darrieus are 
connected to the collar, which is integral with the rotating sleeve, by 
an airfoil shaped piece that slides into the hollow blades and is 
secured with a rivet. See ^figure 13c. for detailed and dimensioned 
drawings. See appendix page D-l for the calculations determining 
appropriate pin sizes. 

GIROMILL 

The giromill has hinged blades that fold up into the center 
shaft when stored. See figure 13a. for stored position. Like the 
Darrieus, upon landing, the lid will open and the giromill will be 
elevated out of i& lander compartment. The giromill struts are 
hinged to the rotating sleeve which allows them to fall into their 
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operational positions once out of the lander compartment. See figure 
12 for elevated working position. A diagonal wire, made of titanium, 
connecting the top strut to the bottom strut keeps the struts in their 
desired positions and prevents them from falling back down into the 
shaft. This wire has a ball and socket connection at the shaft. See 
figure 13d. for detailed drawings and dimensioning. See appendix 
page D-2 for the calculation determining appropriate pin sizes. 

The generator for both of these machines will be connected to 
the stationary shaft by a bracket. 
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Our main goal was to design a functioning vertical-axis wind 
turbine that would be lighter than the existing design of a tornado 
vortex wind turbine. At an estimated shipping cost of $20,000/kg, it 
is obvious why low mass is crucial to the design. The total masses of 
the three machines are listed below. 

Tornado Vortex =30.0 kg ^ 

Darrieus 1.624 kg ) 

Giromill 0.900 kg ^ 


Although the vertical-axis wind turbines we have designed 
would be quite expensive to produce (especially the Darrieus), the 
reduction in shipping cost is so large that the production cost is 
almost negligible. 

We also see that the Giromill appears to be a better choice than 
a Darrieus. The small size of our machines leads to small stresses 
developed in the blades. Because of these small stresses, the ^ 

structural advantage of the Darrieus is not a major benefit. If the 
machine were to be scaled to a larger size, the structural advantages 
of the Darrieus would have a much larger effect. As far as a 1 WatJ 
turbine goes, the Giromill is a better choice than the Darrieus. 

Although the varying cross-sectional Darrieus is not the best ^ 
choice for a 1 Watt Mars turbine, it deserves serious consideration 
for future wind turbines on Earth. By varying the cross-section, both 
the required mass and the maximum stress are significantly reduced. 
Although initial production costs would be high, a mold could be built 
so that several of the turbines could be built at moderate cost. 
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Future Work 
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In the development of this design, several ideas came up which 
we were unable to incorporate because of temporal constraints. 

These included: . , 

•Use of magnetic bearings instead of ball bearings. This would 

cut the mechanical losses considerably, but magnetic 
bearings do take a considerable amount of power. 

Perhaps if a material were created with a critical 
temperature of superconductivity above the range 
which we're working... 

•Use of pitch control, or "smart blades," on the giromill. 

varying the pitch of the giromill blades as they progress 
around the shaft, either by an active control system 
(computer feedback loop and servo motors) or by 
ingeniously locating the pin between the blade and strut 
at a point other than the quarter-chord point, where a 
varying coefficient of moment could be used against a 
spring to change the pitch of the blade. 

•Develop a shock absorber system to mount the wind turbine 
on while stored in the lander. The purpose of this would 
be to cut the acceleration of landing and further reduce 
the size of structural components. 

•Scale up the wind turbine so that it could provide enough 

power to operate the research package in the transmitting 
mode (about eight watts; p. A-2). 

•Wind tunnel test models of the machines to insure that 

overspeed will not occur. While the best information we 
could find has indicated that this will not occur, airfoil 
data does not cover angles of attack over approximately 
20 degrees; in high winds, the angle of attack could range 
from 0 to 180 degrees. If overspeed did occur, we 
recommend the addition of a braking system, either 
aerodynamic or mechanical. Either could be deployed by 
centrifugal force acting on a governor. y 
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Darrieus Machine Configuration 



Giromill Machine Configuration 
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Must operate in storms 
Must operate in Martian environment 
Want a lightweight design 
Want a simple design 
Must withstand landing 
Want sufficient life for mission 
Must produce 1 Watt of power 
Want small design 

Must not interfere with experiments 

Want to be adaptable to changing 
wind velocities 

Must convert energy from wind 

Want to provide power continuously 

Must provide means for starting 

Want dynamic stability 

Want machine to stay upright 

Want to achieve maximum power-to- 
mass ratio 

Must eliminate need for maintenance 

Want to resist corrosion 

Want to use wind from any direction 

Want to minimize friction losses 

Want maximum aerodynamic performance 


Function 

Withstand Storms 
Resist Environment 
Minimize Weight 
Maximize Simplicity 
Withstand Impact 
Optimize Life 
Provide Power 
Minimize Size 
Prevent Interference 
Promote Adaptability 

Convert Energy 
Ensure Continuity 
Ensure Starting 
Ensure Stability 
Maintain Orientation 
Maximize "Power/Mass" 

Eliminate Maintenance 
Resist Corrosion 
Accommodate Orientation 
Minimize Friction 
Maximize Aerodynamics 
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Want to minimize mechanical vibrations 

Want to minimize surface roughness 
( to minimize drag) 


Minimize Vibrations 
Resist Abrasion 
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TABLE 2-3. - Zonal wind cumulative probabilities versus season (m/sec, + from west) 


L, (deg) 
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-12.7 
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- 11.2 

-3.4 

-25 
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9.0 
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11.7 
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2484 
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- Meridional wind cumulative probabilities versus season (m/sec, ♦ from South) 
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A moving mass of air, as shown in the figure below, has kinetic energy of 
motion. For the air moving with the wind velocity, V, the kinetic energy is 

T=0.5mV 2 


where T=kinetic energy 
m=mass of air 
V=velocity of air 



For our parcel of air, the mass is equal to the density times the volume, or: 


m=pAx 


where 


p=mass density of the air 
A=cross-sectional area of the parcel 
x=length of the air parcel 



Substituting this mass into our energy equation yields. 


T=0.5(pAx)V 2 

The power in the air, P a j r , is then given by the derivative of the kinetic 
energy with respect to time. Assuming incompressible flow, we have: 


P a j r = dT/dt = 0.5pAV 2 (dx/dt)+ 0.5pAx(2V)(dV/dt) 
For steady air flow, (dV/dt)=0, and (dx/dt)=V. 
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The power equation now reduces to 


P air =0.5pAV 3 

When wind flows past a windmill or turbine, as in the diagram shown 
below, the maximum power that can possibly be extracted from the wind 
can be calculated. 



Circular tube of air flowing through ideal wind turbine. 


For this idealized case, the air is moving from point 1 to point 4 past 
the turbine. If the air flow is ideal, and the maximum possible power is 
extracted, it can be shown by momentum theory [*] that 

V 2 = V 3 = (2/3)Vj A 2 = A 3 = (3/2)Aj 

V 4 = (1/3)V, A 4 = 3A! 

Now, taking an energy balance, for this idealized case, 

^mech, ideal" * > airl'* > air4 
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where P m ech,ideal = the maximum P ower that can be exlracted 
P a iri = the power in the air at position 1 

P a ir 4 = the power in the air at position 4 
Substituting into our power equation, 

Pmech, ideal = OJp^V ‘ A 4 V 4 3 ) = 0.5p(8/9)Aj V, 3 
Expressing in terms of the physically meaningful terms, A 2 an d Vj , 

Pmech, ideal = 0.5p( 8/9)(2/3 ) A 2 V , 3 = 0.5(16/27)pA 2 V , 3 


where (16/27) is the Betz coefficient. 

The actual mechanical power is typically defined as 

Pmech=Cp (0 - 5pA 2 V l 3) 

where Cp is defined as the coefficient of power and must be less than the 
Betz coefficient. 

Wind speed probability relationships 

The figures shown on the next page show the relationship between 
,X, and the shaft power, where X is the tip speed ratio defined as: 


^ ~ ^ax 0 ^ 


where: 


r mov = maximum horizontal distance 
m ax 

from shaft to blades (*b* in troposkien notation) 
a>= angular velocity of shaft 

= free stream wind velocity 
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curves for other Darrieus turbines all look nearly identical. One 
important thing to notice is that all the curves have a maximum C p of 
about 0.35 at a tip-speed ratio of just less than 6. The other major 
point of interest is that the maximum power generated occurs when ? 
the wind speed is about twice the wind speed for the maximum,^ 
coefficient of power (C pm ). ^ 

In order to determine what speed to operate our turbine at, we 
need to analyze the wind speed probability. Data for wind speed on 
Mars is very difficult to find, and what is found is only for two 
particular locations. Since our wind speed probability function 
depends not only on time, but also on location where we land, the 
global mean wind speed was used. A probability analysis was then 
done using the global wind speed as our base. 

The most important item to determine is the average power we 
are able to generate over time. This can be estimated with a 
probability study. — 
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The previous curves show the relationship between the power 
and the wind speed. At first it seems odd to have the relationship 
between power and wind speed be nearly linear, since power depends 
on wind speed cubed. However, we must remember that efficiency 
goes down with increasing wind speed. For the previous curve, a line 
appears to be a good fit for the curve; however, for other turbines it 
has been foundf-'?] that a better fit can be made with the relationship: 


Pe - 0 u<u c 

P e = a + bu k Uc<u<ur 

Pe - PeR u>ur 


where: 

P e = usable electric power generated 
PeR= the rated (maximum) electrical power 
u » wind speed 

ur = wind speed at which maximum power is generated 
u c =* wind speed when mechanical and electrical losses are 
equal to shaft power 
a,b= constants used to fit curve 
k =» the Weibull shape parameter, a probability term 


It is suggested [6] that if the wind speed probability is not well 
known that a value of k»2 should be used. This is the value we will 
use. The values of a and b that give the best curve fit are 


a = 



u c -u 


k 

R 


b 


_sa 


u R -u 


k 
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With this relationship, the power is now approximated with 
the curve on the following page. 
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Now that we have an equation for power as a function of wind 
speed, we are able to calculate our average power from the 
relationship: 


Pe.avg = J Pef(u)du 
u=0 

where: 

f(u) is the probability density function of wind speeds; it has 

oo 

the property j f(u)du » 1 
u-0 

A recommended probability density function is [£]: 



where c is a constant; c* 1.12umQgn 

Substituting this into the above equation for P e ,avg , yields: 


I - (a 



oo 


Ur 

Pc.avg = J (a+bu k )f(u) du + Pc R j f(u)du 

U C Ur 


This equation can be integrated if the change of variables is made so 
that 



Then, in terms of x, we have 
Jf(u)du = Je _x dx = -e x 


Ju k f(u)du = 



f(u)du 


Jc k xe-*d 


x 


= -c k (x+l )e* x 


Substituting in the limits of integration yields 


Pe.avg - PeR 


K K 

e~(uc/c) e -(uR/c)| 

(u R /c) k - (u c /c) K J 


The quantity in brackets is often referred to as the plant factor (PF) 
The value of u c is nearly always in the range 0.4u R < u c < 0.5u R . 

The normalized power is defined as 


Pn=(PF) 



The following curve shows Pn for u c =0.4ur . Since we chose a value 
of k = 2.0, we see that we have a maximum normalized power if u R /c. 
equals 2.0. By designing our turbine for this situation we will 
produce the most usable energy over time. 
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Since ur/c = 2.0 and c=1.12u me an , we can deduce that ur = 
1.8u me an will be the optimum selection. 

Reviewing the curve of C p vs. X and using the above expressioil 
yields for the optimum case: v ' 

\r*2.9 = 


... O'm/^ 

Solving this equation tells us that the optimum (r m ax(0) = 31.3 m/s. 

With this information, we are now able to to size our Darrieus. We 
decided to get at least 1W of power from wind speeds greater than 6 
m/s. Returning to our power equation, and an assuming overall 
efficiency of rj =0.8 to cover all mechanical and electrical losses, we 
have: 

p v i 

1 Watt = 0.8 { Cp( u =6m/s)(0.5)(0.01665 kg/m 3 )A SW ept(6m/s) 3 } 
where 

Cp( u =6m/s) is taken off the chart for our calculated A.( u= 6m/s)=5.2 1 
We find Cp(u=6m/s) = 0.35. ^ 


Substituting this in yields: 
Aswept = 2.00 m 2 
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List of Symbols 

One-half the total height of our blade 
The area swept out by both blades 
The ratio of the swept area to the arclength 
The maximum horizontal position of our blade 
Centrifugal force 

A constant concerning the varying density 
acceleration of gravity, 3.70 ^ on Mars 

Gravitational force acting on a section of the blade 
the step size used for numerical integration 
a parameter used in integration 
a variable used in bisection method to solve for K 
a variable used in bisection method to solve for K 
tip-speed ratio 

an arbitrary point used in the derivation, see 
Figure 1 

Horizontal coordinate 

The average horizontal coordinate for the blade 

The mass density per unit length of our blade 

The mass per unit length of a thin airfoil "skin" 

The portion of the density (per unit length) 
that is constant 

The density (per unit length) of the internal blade 
support @ z=0 

The portion of the density (per unit length) that is 
varied 

Length of blade between point of maximum 
horizontal deflection and point P, see Figure 1 

Total arclength of one blade 
Tension at arbitrary point P 
Tension at z=0 

Slope of blade at point P, see Figure 1 
The angular velocity of our blade 
Rotational parameters, see Eq. (11) & (12) of 
derivation 
A ratio of densities 
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Referring to Figure 1, we obtain two equations that must be 
satisfied for equilibrium. 


IF Z = 0 


IF r = 0 


where IF Z = sum of all forces in z direction 

LF r = sum of all forces in r direction 

For our situation these equations reduce to 

( 1 ) TsinG = Cf ^ 

(2) Tcos9 = T 0 + G 

where T = tension in member 
0 = angle in Figure 1 
T 0 = tension at vertical midpoint 

s 

G = gravity force = Jpg d s 

0 

s 

Cf = centrifugal force = J p<o 2 r d s 

0 

where p = mass per unit length 

0) = angular velocity 
s = arc length 
g = acceleration of gravity 


dr 

Taking the ratio of equation (1) and (2) and noting that tan0= - 


Cf dr 
tan 0 = T 0 + G - ' dz 


Substituting in to equation 3 yields 


T- »Z 


r 



(4) 


dr _ 
dz ~ 


s 


J pco 2 r ds 

0 


T 0 + G 


Equation (4) is subject to the boundary conditions 

r = 0 at z = a ^ 
dr 

~r~ - 0 at z = 0 ✓ 

dz 

Assuming a rotational speed of about 40 rad/s and considering any 
point with radial position of greater than 0.1 meters, (our blade's 
average radius is 0.677 meters) 


Centrifugal acceleration = co 2 r > (40rad/s) 2 (0.1m) =160 


m 

sec 2 


Gravitational acceleration = 3.70 


m 

sec 2 




Clearly, the gravitational acceleration can be neglected. After we 
select our material, we will also show that the aerodynamic forces 
are negligible in determining the shape. 


For constant rotational speed, equation (4) reduces to 


(5) 


dr co 2 r 
d~z = -T 0 J^ Td 


0 




When I observed a conventional troposkein solution from Sandia[l], ^ 
noticed that the tension varied according to the equation 


( 6 ) 


T r 2 

T 0 = 1 * C ( b 2 _1) 





where C was some constant 


In observing equation (6), it was observed- that since the tension 
varied along the length of the blade, it might be a good idea to vary 
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some of our mass density so that we may reduce mass where it is not 
needed. This will not only reduce the mass in that location, but will 
lower the stresses throughout the blade. Remember that density is 
density per unit length, so what is really being varied is the cross- 
sectional area. / 

? 

Therefore, it was decided £s thougly we might have some portion of 
our blade cross section wiuT'T'cbnstant cross-section(i.e. a thin airfoil 
"skin"), and then an internal support with a varying cross-section. 

Remember that density is defined as mass per unit length, so a 

varying density really means a varying cross-s ection al— area. , » * /, 

The total density would then be in the form > 

rl / 

P = Paf + Psupp(l'D(^ *1) ) 

where 


Paf = the density of the airfoil skin; a constant 
Psupp = the density of the internal support afo z=0 
D = a constant to be optimized 


6 ^ 




/ 


Eventually we decided to make our airfoil as one unit, and for our 
case paf is equal to zero. The derivation is for the more general case 
and may be applied to our case by just setting p a f=0. 


The density equation may be rewritten as 

xl 

P = Paf + Psuppd+D) - Psupp(Djj2 ) 


For simplicity these will be grouped in to two terms, one which does 
not depend on r, and one that does. The density terms will be: 

p c = the constant density portion 


p v = the varying density portion 
Our density now will be written as 
(7) P = Pc - Pvt 2 


T- 1 + 


r 



Putting equation (7) into equation (5) yields 


( 8 ) 


dr 

dz 


co 2 p c 


J rds 


0) 2 p v 


T 0 0 J‘~ ' T 0 


r 3 d s 


0 


We can rewrite equation (8) by noticing that 
(9) ds 




) 2 d z 


( 10 ) 


where 


dr _ . 

3“ = -^c 2 
dz c 


z ^ 

M 1+ (^) 2dz +Q v 2 Jr 3 "\| l+(^) 2 d 

0 0 


(11) 

CO 2 Or 
2 _ UJ HC 

*0 

(12) 

. C0 2 Pv 

Q 2 — r 

*0 



Now, we can change the integro-differential equation (10) to an ordinary 
differential equation 


(13) 


d 2 r 

dz 2 


= 1+ffi + a,V-^ l+(^)2 


After some algebraic manipulation, equation (13) can be written in the 
form of exact differentials. 


(14) 


d_ 

dz 


V i+ & 2 =- 


Ov 2 d 


dr.. «£l_ fr 2 x + ifvi ^ 

2 dz (r ' 4 dz 


dr 


/ 


Integrating and substituting in the boundary condition that • — 0 at r=b 


T-ir 


r 



(15) -\j 1 +<£)2 = 1 - (r 2 -b 2 ) + ^ (H-b“) 

Squaring both sides yields 

(16) = [1 - ^ (r 2 -b 2 ) + ^ (r*-b«) ] 2 


Evaluating the left hand side (L.H.S.) yields 

tuc _ 1 , / Q c 2 b 2 zL Qy 2 b 4 T 4 

L.H.S. - 1 + ( 2 V _1) ' 4 ( b 4 _1) )2 


Qy^b 4 r 4 

-A 2 ( b 2 -1) - 4 


<M -1) > 


This may be rewritten into equation (16) as 

( dj\ ? _ ( flc 2 b 2 sL n v 2 b 4 r 4 X/OoVrl fi„ 2 b 4 T* \ 

W ( 2 V* -•> • — X "V" V '> - ~V" £ -‘) -i) 

This equation may now be rewritten as 


^.d-T ^ Oy*b g / l£s r 2 ti. \ / ^Pc r2 r 4 o X 

^dz^ - 16 ''Pvb 2 ^b 2 ’ 1 ^ * ^b 4 Vp v b 2 ^b 2 ' 1 ^ * ^,4 *!) ' Q b 4 ) 


dr 2 _ n/b» _ s /ik_ r 2 *1 \ / 

xlz' 16 n v 2 b 4 ^Pvb 2 ( b 2 ‘ 1} ' <b 4 •*>) ( 


2Es_ r 2_ 

Pvb 2 < b 2_1) ‘ ( 


b 4 


1 ) 


Q,b 4 


The above equation was obtained by noting that 



Ac 2 


Qy 2 


£c 

Pv 


Now we will define two new parameters 



?c " O J ' & 


n v 2b 4 

8 


K2 


T- IL> 


Now we may write 


(£) 2 = 4K2 (£(£-» • V -1>)(K2$£-1> - £ '!>) -0 
Taking the square root of both sides yields 


^ = ±2K*\J( ^-D - <£ -d)(k 2 (^-d - -n) - 1 ) 


. dr 

If we observe Figure 1, we can see that £ j z 


is going to be negative. 


j£ = -2K-^( £<£•!) • •1 ))(k2(^-D -<£ •'» -0 


(17) 


We may now solve equation (17) numerically with its two boundary 
conditions for both the shape r(z) and the unknown K 


?* 


T- 1 7- 
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***** ********* **** This prograa will generate solutions for a varying density troposkien ************** 
PROGRAM TROPO SKEIN 

****************** T00MIGH and TOOLOW will be used to determine the correct value for K ************** 


LOGICAL TOOHIGH, TOOLOW 



A is one-half the height of the blade 
2 and R contain the coordinates of the blade 
H is a step sire used for numerical integration 
K is what I called "K squared" in the derivation 
KLO and KHI will be used to solve for K using a bisection Method ***** * * * ** *** ****** 



REAL RPRIME, A, 2(1000), R(1000),H,CRCY<1000) 
ft , CRCX(1000) / K ( ,KLO / KHI / HAXTEN / TENSARRAY(IOOO) 
ft , RATI0(31 ) , ADIVBC31 ), 2DIVAC1000) , STRESSARRAY(IOOO) 
REAL DERIVATIVE (1000) , FI, F2, F3, F4,HTIMESD 
EXTERNAL RPRIHE, AREASVEPT, RAVG, ARCLENGTH 


****************** Write the output to a file * ******* *** * A k A ****************************** * * kkk *** *** 

OPEN (18,f i le*' tropoutl 1 ) 

UR1TE(18,55) 

WRITE(18,*> 

WRITEC18,*) 


This nain loop will carry out the entire calculations for varying values of A and B ** 
so that the area swept is nearly constant. 


LOOPSTEPS * 31 

C DO 1000 LOOPER = 1,L00P$TEP$ 

C A 3 0.55+ (REAL(L00PER-1))*0. 35/30.0 


***************** Set B so the area swept is nearly a constant for coaparison's sake ****************** 
C 1.31/(2*A> 


If we want to work with a particular a and b, we can use the following 
****************** two lines and coaaent out the "DO 1000 loop" 

A* 0.737 
B=* 0.889 


*********** Our boundary condition requires the next two lines to be true *********************** 

Z(1) 3 A 

R(1)*0.0 


*********** Next I set a step sire and a "guess" for the value of K 
H= A/800.0 
K*1 .4 
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******* ** ** Set KLO and KHI to the upper and lower bounds for K *** * * ** *********************^ h** 

KLO-0.0 

KH 1*45.0 


NPTS*800 


NPTS is the nuaber of points I will use to approxiaate the integral **************** 


The next line is required so that our boundary condition is aet. 


DERIVATIVE(NPTS)*0. 0 


HAXSTEPS is the aaxiaua M of tiaes I will carry out the bisection aethod *********** 


HAXSTEPS * 30 


100 


This following loop will carry out the bisection aethod to solve for K ************ 

DO 33 HMM * 1, HAXSTEPS 
T00HIGH*. FALSE. 

TOOLOW -.FALSE. 

This loop will carry out the integration using the fourth order Runge-Kutta ********* 
Method and keeps the values of r and 2 in the arrays R and Z *********************** 
DO 100 1*2, NPTS 

Z(I) * A - (REAL(I))*A/REAL(NPTS) 

F1*H*RPRIME(R(I-1),K,B) 

F2*H*RPRIME(RC 1-1 )+0. 5*F1 , K, B) 

F3=H*RPRIHE(R(I-1)*0.5*F2,K,B) 

F4=H*RPRIHE(R(1-1)+F3,K,B) 

HTIHESD is the step size tiaes the derivative **************************************** 
HTIMESD * (F1 + (2.0*(F2+F3))«’F4)/6.0 
R ( I ) * RU-1) + HTIMESD 
DERIVATIVE (1-1) - HTIMESD/H 
CONTINUE 


**************** End of integration loop 


If our slope at the end does not aatch our boundary condition, aake K higher ****** 
IF (DERIVATIVE (NPTS-1 ) .GT. 0.0001 ) THEN 
KLOK 

K * (K+KHD/2.0 
TOOLOW * .TRUE. 

END IF 

*********** gnd of slope at boundary equation check 


If we select too high a value for K, we will get an iaaginary derivative, which *** 
is not a solution to our p rob lea. The follwing check prevents these solutions 


**************** froa appearing, and corrects K so that the next atteapt will be better. 
IF (.NOT. (TOOLOW)) THEN 

IF (DERIVATIVE (NPTS-2) . EQ. 0 . 0) THEN 
KHI*K 

K*(K*KL0)/2.0 
TOOHIGH » .TRUE. 

ENOIF 

ENPIF 


f 

s 

h 

\ 

i 
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** ************** End of the check for K too high *************************************************** 


r*********** if our K Is now the solution, stop performing the bisection method 
IF ((.N0T.(T00HIGH)).AND. (.NOT. 

I (TOOLOW))) THEN 
GOTO 444 


END IF 

**************** 


End of check for correct K 


** ** * **** ******* 
33 CONTINUE 

* 1 ■ * I l lLLllll lj L 

A"e«A XXAee AwAAA A 

444 CONTINUE 


If K Is incorrect, go back with an updated value for K * ** * * * * * * * *** *** ** * ** *** * * * 


*************** now that we have the solution, calculate the desired parameters ******************* 
*************** and print out the output ********************************************************** 

ASPT * AREASUEPT ( R, H,NPTS) 

ARCLT * ARCLENGTH(R, DERIVATIVE, NPTS,H,K,B) 

HAXTEN * $QRT< 1.0 + (DERIVATIVEC1 ))**2) 

URITE(18,999) A^^ASPT^RCL^RAVGCR^PTSMAm^ASPT/ARCLT 
RATIO(LOOPER) * ASPT/ARCLT 
ADIVB(LOOPER) * A/B 
1000 CONTINUE 


End of main loop for different values of a and b 


For a given a and b, the following loop 
array with the values of the tension at 


DO 93 LFT * 1,NPTS 

TENSARRAY(LFT) * SQRTd.O + (DERIVATIVE(LFT)>**2> 
ZDIVA(LFT) * Z(LFT)/A 
93 CONTINUE 


will produce an 
every point 


End of tension calculating loop 


This loop gives the tension and 
us the total mass required with 


stress ratios at any given point as well as giving *** 


RHOSUW-O.O 
DO 94 LFT2 ■ 1,NPTS 
VARY-1. Of1.55*C1-<R(LFT2)7B)**2) 
STRES$ARRAY(LFT2)*TENSARRAY(LFT2 ) /vary 
RH0SUN-RH0SUN+VARY*(1*SQRT(1+ ( DERIVATIVE(LFT2) )**2) )*H 
94 CONTINUE 
C VRITE<*,*> RHOSUN 



HNPTS * 51 


The following loop will enable us to plot a circle of the same 
arclength as our troposkien for the purpose of comparison *** 


T-2-0 
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PI * ACOSl-1.0) 

CRCY(HMPTS) « ARCLT/PI 
CRCX(NHPTS) * 0.0 
DO 9432 LCC » 1,HNPTS-1 

PHI » REAL(LCC-1)*PI/(2.0*<REALO«1PTS))) 

CRCY(LCC) « (ARCLT*SIN(PHI))/PI 
CRCXCLCC) * (ARCLT*COS(PHI))/PI 
9432 CONTINUE 

End of loop to calculat* clrela coordinates 


*************** 


End of the entire calculation process 


*************** These lines can be used to plot various parameters *i 
C CALL XYUNITl ' meters' ,' meters' ) 

C CALL CURVICRCX, CRCY,MHPTS, ' CIRCLE ' , ' ' ,1, .TRUE. ) 

C CALL CURVIR, Z,NPTS, 1 TROPOSKIEN' , 1 ',3,. FALSE.) 

C CALL CURV(ADIVB,RATIO / LOOPSTEPS / ' Asuept/S' , 1 ',2,. FALSE.) 

C CALL SPLOTCratio.ps' , 'POUER-TO-MASS' , 'a/b' , 

C £ ' AREA/LENGTH 1 , ' t ' ,5 .0,5.0, .TRUE. , .TRUE. ,4,4, 

C & 5,5,0.6,1.5,0.69,0.75) 

C CALL SPL0TCtskein2.ps', 'COMPARISON', 'R', 

C £ 'Z', 'f , 5.0,5. 0,. TRUE.,. TRUE. ,4,4,5,5, 

C £ 0.0,1. 0,0.0,1.0) 

C CALL CURV ( ZDIVA, TENSARRAY, NPTS, 1 TENS ION RATIO' ," ,2, . FALSE.) 

C CALL CURV(ZDIVA,STRESSARRAY,NPTS, ' STRESS RATIO' ," ,3, .FALSE.) 

C CALL SPLOTC tens. data', 'TENSION AND STRESS RATIOS' , 'Z/A', 

C £ 'RATIOS', 't', 5. 0,5.0,. TRUE.,. TRUE., 5,10,7,10,0.0,1.0,0.5,4.0) 

55 FORHAT(//,9X, ' A ' , 12X, 'B* ,12X, ' K' , 11X, ' Aswept' ,10X, 'S' , 

£ 13X, 'R',7X, ' (T/To)max‘ ,6X, 'Asuept/S' ) 

999 FORHAT(5X,E9.3,4X,E8.3,4X,E10.4,4X,E10.4,4X,E10.4,4X,E9.3, 

£ 4X,E10.4,4X,E10.4) 


STOP 

END 


END OF THE MAIN PROGRAM 


Function RPRIHE calculates dr/dz at a given point 
****************** The variable X in thi a function is what I call r in the rest of the program ***** 
FUNCTION RPRINE(X,K,B) ' 

REAL X,B,Q,L,H,K, ZETA, PRODUCT, ROOT 


****************** Zeta will be changed manually to optimize 
ZETA « 1.3 


****************** The following are parameters to save calculation time 
Q ■ 2.0*ZETA/(B**2) 

L - «X/B)**4)-1.0d0 
H » (<X/B)**2) - I.OdO 
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PRODUCT * (Q*N) - L 

ROOT * PRODUCT*((K*PRODUCT)-T.O) 


** *************** * 


If K is too high, root will be negative, and this will lead to imaginary 
solutions, which nake no sense for this probtea. 


IF (ROOT. GE. 0.0) THEN 

RPRIHE > SQRTCK/4.0) * SORT (ROOT) 


***************** If K it too high then ROOT will be negative. Too avoid coaputer error, set ***** 
***************** RPRIHE * 0. Then in the aain prograa if RPRIHE * 0, I change K so that it **** 
***************** wit l be lower for the next Iteration. 

ELSE 

RPRIHE - 0.0 


END IF 

RETURN 

END 

******************************** END OF FUNCTION RPRIHE ********************************************** 


**************** Function AREASWEPT will calculate the area swept by both blades ********************** 
FUNCTION AREASWEPT ( R, H, NPTS) 

REAL R(NPT$),H 
AREASUN >0.0 
DO 5 III=1,NPTS 
AREASUH > AREASUH +R(III)*H 
5 CONTINUE 

AREASWEPT > 4.0*AREA$UH 

RETURN 

END 


END OF FUNCTION AREASWEPT 


Function RAVG will calculate the aean value of r for our blade ** * ** *** * ***** ****** *** 


FUNCTION RAVG (R, NPTS) 

REAL R(NPTS) 

SUMAVG >0.0 
DO 6 JJJ>1,NPTS 
SUHA VG> SUHA VG * R(JJJ> ' 

6 CONTINUE 

RAVG * SUHAVG/ REAL (NPTS) 

RETURN 

END 

**************************************** END OF FUNCTION RAVG 


T- 


f 
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****************** Function ARC LENGTH will calculate the arclength required of each blade * *** * **** *** 
FUNCTION ARCLENGTH(R, DERIVATIVE, NPTS, H, K, B) 

REAL ARCSUM, R(NPTS) , DERIVATIVE (NPTS) , H, K 
ARCSUH * 0.0 
DO 7 LLL“1 , NPTS 

ARCSUN * ARCSUH+(SQRT(1+(DERIVATIVE(LLL) ) 

I **2))*H 

7 CONTINUE 


ARCLENGTH * 2.0*ARCSUH 

RETURN 

END 



I 

!■ 

* _ 

l 

r_ 

r 

r_ 
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Discussion on computer sol ution to troooskien equations 


We begin with the problem of trying to solve equation (17) 
together with its two boundary equations. Since it is a first order 
differential equation with two boundary conditions and one 
additional unknown, we see that we have the right number of 
equations in order to solve for both the shape ( r(z) ) and the 

unknown K. The Runge-Kutta method of order 4 and the bisection 

method are applied simultaneously to equation (17). The procedure 
is outlined below. 

First, several values of a and b are computed by the computer 

so that all of these values sweep out roughly 2 m 2 (the required size 

for 1 Watt of power). The solution is then generated for all of these 
values of a and b so that the optimal height-to-width ratio will be 

obtained. Also a value of £ is selected and this value will be changed 

manually after observation of results in order to optimize the value 
of this parameter. Now, for each combination of a,b, and z the only 
unknown left in equation (17) is K. ^ 

First we make a guess value for K, and set upper and lower 

bounds for K. The lower bound is 0. The upper bound and the guess 

value are found after trying the solution a couple of times. Then the 
equation is numerically integrated with this value of K. 

The numerical integration is done with the fourth order Runge- 
Kutta method. The method can be found in any numerical 
mathematics textbook and may be seen clearly in the included 
computer program. The step size was continually decreased until no 
change in the results resulted from a further decrease in the step 

size. The required step size was found to be . 

In order to solve for K, we must ensure that our solution 
satisfies the boundary conditions of no slope at z=0 and r=b at z=0. If 
K is selected too low, we run into two problems. First, the solution 
does not reach b when z=0. Second, there is a finite (non-zero) slope 
at z=0. Clearly, too low a value for K is not a solution to the equation. 
Also, if K is selected too high, the solution "steps" past b and the 
product under the radical in equation (17) will be negative. This 
violates our second boundary condition and is also not a solution. 

This brings us to the bisection method used. 

After attempting the calculations with the guess value for K, 
we could determine whether K was too low or too high from the 
results. After making this determination, K was updated using a 


r ' 


bisection method. The method is outlined below. 

Initial Guess of 
value for K 


KLO = 0 = Preset 
lower bound for 
K 


Actual Value 
of K 


Step 1 


KHI = a preselected 
value for an upper bound 
for K, determined by trial- 
and-error 


We found that K was too high, so update its value 
and the value of KHI 

I 


— 

| 

New 
guess 
for K 

i 

i 

i 

Updated value for 

KH, | 1 

— 

1 

KLO 


Actual 





Value 


F 1 



of K 



Step 2 


Now, we find that K is too low, so update K and KLO 


New 


New 

| guess 


1 

■ 

1 

1 

1 

[ 1 

Actual 
Value 
of K 

1 

1 

KHI 


Step 3 


The procedure continues and usually converges after about 8 
steps, since our initial interval is quite small. 


f- 


i 


~r~ z<r 


The whole procedure is carried out for 30 different values of a 
and b for each selected Different solutions are then obtained for 
different values of £. 

The area swept is then calculated by the formula : 

a 

A swept = 4 J r d z 
* o 

This is done numerically by noting that 

NPTS 

A swept “ 4 Az J] r i 
i=l 


where NPTS = T — 
A z 


The arclength is calculated from the formula 


= 2 W 


1 + (^) 2 dz 


This is solved numerically with the approximation 
NPTS 

I 77 

\2 


S - 2AzX V I + 

i*l 


Also calculated in the computer program is the average 
horizontal location along the blade; it is called RAVG or R. 

The most important parameter generated is A swc pt/S, which is 
the ratio of the area swept to the arclength. It is effectively a 
power-to-mass parameter; optimization, of this parameter is vital. 




r 


Everything has now been solved for except for the tension. 
This is done by recalling equations (1) and (2) from the derivation. 
If these two equations are squared and added together, we obtain: 

T 2 = Cf 2 + (T 0 +G) 2 

Dividing all terms by (T 0 +G) 2 we obtain 

T 2 _ -C£- +1 

(T 0 +G)2 ' (T 0 +G)2 + 


Also, from equation (3) in the derivation, 

..H, - (.— )2 + 1 

(T„+G)2 1 dz ’ 


,Zl^' 


Now, gravity can be neglected in the equation, and we obtain: 


T 

Ti 




This can be evaluated numerically at every point. 


Reading Tabular Data 


The following is a brief outline of what all the terms on the tabular 
data mean and how to interpret the results 



A one-half the height of the blades 

B the maximum horizontal position of the blades 

a rotational parameter, no obvious physical 
significance 

A swept the area swept by both blades, power generated is 
directly proportional to A swcp i ^ L ’ 

S the arclength of one blade, effectively a measure 

of required mass 

R the average horizontal coordinate of our blade 

(T/T 0 )max the maximum normalized tension (T 0 is where 
tension is minimum) 


A swept/ S the most important column in the tabular data; a 

power-to-mass ratio; maximizing this parameter 
is crucial to our design 


~T-2 ? 


r 
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'I-- 1,3 


Aswept 


(T/To)max 


0.550E+00 

0.562E+00 

0.573E+00 

0.585E+00 

0.597E+00 

0.608E+00 

0.620E+00 

0.632E+00 

0.643E+00 

0.655E+00 

0.667E+00 

0.678E+00 

0.690E+00 

0.702E+00 

0.713E+00 

0.725E+00 

0.737E+00 

0.748E+00 

0.760E+00 

0.772E+00 

0.783E+00 

0.795E+00 

0.807E+00 

0.818E+00 

0.830E+00 

0.842E+00 

0.853E+00 

0.865E+00 

0.877E+00 

0.888E+00 

0.900E+00 


.119E+01 

. 117E+01 

.114E+01 

. 112E+01 

.110E+01 

•108E+01 

. 106E+01 

.104E+01 

■102E+01 

■100E+01 

-982E+00 

■966E+00 

•949E+00 

-933E+00 

■918E+00 

•903E+00 

.889E+00 

■875E+00 

.862E+00 

•849E+00 

•836E+00 

•824E+00 

■812E+00 

-800E+00 

. 789E+00 

.778E+00 

.768E+00 

.757E+00 

-747E+00 

.737E+00 

.728E+00 


0.3138E+02 

0.3904E+02 

0.4500E+02 

0.3521E+02 

0.2256E+02 

0. 1666E+02 

0.1311E+02 

0. 1068E+02 

0.8915E+01 

0.7574E+01 

0.6509E+01 

0.5658E+01 

0.4955E+01 

0.4375E+01 

0.3885E+01 

0.3465E+01 

0.3108E+01 

0.2800E+01 

0.2528E+01 

0.2291E+01 

0.2081E+01 

0.1898E+01 

0.1735E+01 

0.1589E+01 

0.1459E+01 

0.1343E+01 

0.1236E+01 

0.1142E+01 

0.1055E+01 

0.9775E+00 

0.9064E+00 


0.2054E+01 

0.2215E+01 

0.2336E+01 

0.2339E+01 

0.2266E+01 

0.2218E+01 

0.2181E+01 

0.2151E+01 

0.2127E+01 

0.2106E+01 

0.2087E+01 

0.2070E+01 

0.2055E+01 

0.2042E+01 

0.2030E+01 

0.2018E+01 

0.2007E+01 

0.1998E+01 

0.1988E+01 

0. 1979E+01 

0. 1970E+01 

0.1962E+01 

0.1955E+01 

0. 1948E+01 

0.1941E+01 

0.1935E+01 

0. 1928E+01 

0.1922E+01 

0.1916E+01 

0.1911E+01 

0.1905E+01 


0.2795E+01 

0.2915E+01 

0.2999E+01 

0.2979E+01 

0.2896E+01 

0.2838E+01 

0.2791E+01 

0.2753E+01 

0.2719E+01 

0.2690E+01 

0.2663E+01 

0.2639E+01 

0.2618E+01 

0.2599E+01 

0.2582E+01 

0.2566E+01 

0.2552E+01 

0.2540E+01 

0.2529E+01 

0.2519E+01 

0.251 1E+01 

0.2504E+01 

0.2498E+01 

0.2493E+01 

0.2489E+01 

0.2486E+01 

0.2483E+01 

0.2483E+01 
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Evaluation of Parameters from troposkien solution 


We have solved for the troposkien shape in terms of the 
parameter K. However, this parameter gives us little insight as far as 
the magnitudes of the stresses. In order to obtain these figures we 
must go back and determine the physically meaningful terms in 
terms of K. Recall that 


K 2 = 

Ov 2 


fly 2 b 4 

= 8 

_ C0 2 pv 
' To 

_ia ggyfel 
lo " 8K 2 


( 1 ) 


We also know that 

bco = 31.3 m/s 

(from our tip-speed ratio) 


p v = (by definition) 


Substituting these into Equation (1) yields 
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We are able to obtain an expression for D also, 


c = 



b 2 Psupp 


D 


(3) 


This relationship is also by definition (see derivation) 
A third relationship is that for the maximum tension 


T a - T max — To(T/To)niax 


The optimal case would be if the stress at all points was the 
same. Therefore, we will try to find a relationship between the 
tension and the stress. T m ax is also referred to as T a because it 
occurs when z=a. 

The tensile stress is simply the tension divided by the cross- 
sectional area. Since our density terms are density per unit length, 
they are really the mass density of the material times the cross- 
sectional area. 



T a /A a 

T 0 /A 0 


^T/T 0 )max 

Pa 


where: 


<r-~' 


£ 




a = stress 

A = cross-sectional area 
U means "is proportional to" 

Since our airfoil has no separate "skin", p a f=0 (see derivation) 
and we may replace ”(d" by "=". 


— ~ . p. (T/T 0 )max (4) 

Gq 1+L> 

Immediately we see that the stress ratio is always less than the 
ratio of the tensions which is what was desired. 

We now have 4 equations to work with which enable us to 
solve for the tensile forces, the constant D, and the densities per unit 
length. The equation for £ has been simplified for the case where 
p a f=0. They are presented below as a summary 

122.46p SU ppD 

( 1 ) lo = K 2 

(2) C = b2(l+5> 
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(3) 


Ta — Tjnax — To(T/To)max 


(4) ^ = !+ D (T ' To)m “ 

We see that equation (2) can be solved for D once £ has been 
picked. In order to determine £, the computer program was run 
several times for various choices of £, until we could obtain a ratio in 
equation (4) as close to unity as feasible. If £ was selected too low, 
the tension ratio became very large and the ratio would not approach 
1. If £ was selected very high, we lose the benefits of varying the 
cross-sectional area and can not approach a ratio of 1. Therefore, 
when trying to obtain a ratio in equation (4) close to unity while also 
minimizing arclength, the optimal value of £ was found to be £=1.3. 
There is no closed form solution to show that this is the optimum, 
however iterative computer solutions reveal that this is the case. 

With £=1.3, the computer program was run, and the output of 
interest here was: 


b= 0.889 meters 
(T/T 0 ) ra ax = 3.928 
K2 = 3.108 

With these values known, we are able to solve for D and our 
stress ratio. 


D = 1.55 


We see that we have taken a loading ratio of nearly 4 and 
reduced it to a stress ratio of approximately 1.5. This enables us to 
reduce our support mass significantly. 

Next, we would like to approximate how much better our 
solution is than the traditional constant cross-sectional area 
troposkien solution. To do this, we estimate the mass required for 
the blades of both machines. Remember that our solution provides 
benefits in three separate ways. First, we reduce weight by varying 
the cross-section along the blades. Second, we obtain a better shape 
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than the constant cross-section solution, therefore requiring less 
weight. Finally, since the stresses are mass dependent, the mass we 
removed where it was not needed helps to reduce stresses 
throughout the blade. Since there is no benefit in completely 
developing a constant cross-section solution, we will consider each of 
the three parts separately, and add their contributions to 
approximate our total mass and stress reduction. 

First, we generated a computer solution for the case of £=1-3 and for 
the case C=20 (nearly constant cross-section). The arclength for the 
z=20 solution was longer than the z=1.3 case because of the shape 
benefits of varying the cross-section. The constant cross-section 
solution would need to have a large enough cross-section to support 
the maximum tensile load. The maximum tensile load would also be 
greater for the constant cross-section solution due to more mass 
present, but this topic will be discussed in section (II). Therefore its 
cross-section will be approximated with the cross-section of the 
varying cross-section solution at r=0. 

p(r) = psupp(l+D) constant cross-section solution 

p(r) = Psuppd-D( ^ - D) varying cross-section solution 

To find the total mass of each blade, we integrate the density times 
the differential arclength. 

s 

Mass = J p(r)ds 
o 


We know that ds = y 1 + ( ^p 2 

Now the mass equation can be integrated numerically from z=0 to 
z=a, and the ratio of the constant cross-section to the varying cross 
section can be obtained. This was done and the results are 

MASSconstan t cross-section _ ^ 

MASS varying cross-section 
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We see that varying the cross-section allows us to reduce mass by 
34.2%. 


Step (II) 

The maximum tensile loads were compared for the constant cross- 
section and varying solutions and we found that the constant cross- 
section troposkein had a 54% greater maximum tensile load. Since in 
step (I), we assumed that the two solutions had the same cross- 
section at the point of maximum tension, we see that the stress is 
54% greater in the constant cross-section solution with this 
assumption. 


Conclusion 


Therefore, we find that by varying the cross-section we are able to 
reduce blade mass by over 34%. We also have a maximum stress 
54% less than that of a constant cross-section solution. We conclude 
that varying the cross-section has some remarkable advantages and 
should be implemented in our design. ^ ^ 
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TENSION AND STRESS RATIOS 


Scale: 

H. inch = 0.2000 
V. inch = 0.7000 


Legend: 
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FIND MINIMUM DIAMETER OF INSIDE SHAFT FROM BENDING STRESS 


Fa 24- N 
1 a 0.75m 


For AlB: <nilt a l.i-GPa 


Dinside 


'4- F I 1 


v- oult 


1 

3 

•2 


Dinside - 5.503- ntn 

Hf® a large FS, say 1.75, because the ultimate tensile stress 

SknoSS*?' hS yi6ld strength wil1 be somewhat lower than this, but is 

Dinside- 1.75 - 9.631- mn < — use 10-itm shaft 

*™*********************************************^^f*A 5 * 6 A********* 

OUTER ROTATING SHAFT ************** 

Inner radius of the rotating shaft is dependent on the size of the 

otIS 9 be f W Sch t h ^ a the f^ry shift, we ^ rKa?i„nrcm 
oiLF (p. ) which has an outside diameter of 19nm. 

Mass of two blades: MassBlades a 2- kg 

Acceleration during landing inpact: Acc a 7.9.91- ^ ^ 

Inner radius of shaft: Ri a 9.5-mn l sec 


BENDING: 


Ro 


MassBlades - Acc 
> <rult-T 


2 

+ Ri 


Ro =* 9.502- mm 


BUCKLING: 


2- Ri 

thickness : » 

95.5 

Ro : - Ri + thickness 
Ro = 9 . 699 • nm 


So use an outer diameter of 19.5 nm based on critical buckling stress. 


S - 10 


r 


The Student Edition of MathCAD 2.0 


For Educational Use Only 


m = 1L kg a 1M sec a IT 

mm = 0.001- m 

giromill 

Radius of Giromill : 

Length of Giromill Blades: 

Mass of one strut frcm -R to +R: 
Mass of two blades: 

2 

2- MassStrut- h 

Jgiro : = 

12 

2 

Jgiro * 0.233- kg- m 


kg- m N 9 

Pa a — GPa a io 

2 2 
sec m 


R := 0.55-m 
h :=» 1.0-m 

MassStrut : = 0.020- kg 
MassBlades 0.380- kg 


2 

+ 2 ■ MassBlades - R 


DARRIEUS [approximate blade 

Radius of Darrieus: 

Height of Darrieus: 

Mass of two blades: 

Inner radius of blades: 
Outer radius of blades: 


shape as circular] 

R :« 0.75-m 

h 2-R h - 1.5- m 

MassBlades :« 0.989- kg 
rl :- 0.75-m 
r2 :- 0.755- m 


Jdarr 


2 2 
rl + r2 

MassBlades- 


4 


2 

Jdarr =* 0.28- kg- m 


S-n 

r 
i ■ 


• Pa 


6ingle Row Deep Groove Ball Bearings 


PLAIN 
— B — 


TYPE Z 


TYPE 2Z 


TYPE RSI 


TYPE 2RS1 


o 

0 - 

a 


a 

^■i 


78? 

■ 


L 

o 


1 

o 


pj 



i 

gp 



ONE SHIELD 


TWO SHIELDS 
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Principal dimensions 


Basic 

load 

rating 

dynamic 

Speed 

rating 1 

Max. 

fillet 

radius 

Mass 



Designation 





d 

D 

B 

C 




Open 

One shield 

Two shields 

One seal 

Two seals 

mm 

In. 

mm 

in. 

mm 

in. 

N 

Ibf. 

rpm 

mm 

in. 

kg 

lbs. 






7 

0.2756 

14 

0.5512 

5 

0 1969 

753 

169 

63 000 

0.15 
0 006 

0.0020 

0.0044 

- 

- 

AX7ZZ 

- 

- 


17 

0.6693 

5 

0 1969 

1170 

263 

56 000 

0.3 

0012 

0.0049 

0.011 

AY 7 

- 

AY7ZZ 

- 

- 


19 

0.7480 

6 

0.2362 

1 720 
387 

40 000 

0.3 

0.012 

0.0075 

0.016 

607 

607 Z 

607 2Z 

607 RSI 

607 2RS1 


22 

0.8661 

7 

0.2756 

3 250 
731 

36 000 

0.3 

0.012 

0.013 

0.028 

627 

627 Z 

627 2Z 

627 RSI 

627 2RS1 

6 

0.3150 

16 

0 6299 

4* 

0.1575 

1 040 

234 

56 000 

0.2 

0.008 

0.0031 

0.0068 

X8 

- 

X8ZZ 

- 

- 


19 

0.7400 

6 

0.2362 

1 460 
328 

50 000 

0.3 

0.012 

0 0071 
0.016 

AYS 

- 

AY8ZZ 

- 

- 


22 

0.8661 

7 

0.2756 

3 250 
731 

38 000 

0-3 

0.012 

0.012 

0.026 

608 

608 Z 

608 2Z 

608 RSI 

606 2 RSI 

9 

0.3543 

17 

0.6693 

4* 

0.1575 

1 110 
250 

53 000 

0.2 

0.008 

0.0034 

0.0075 

X9 

- 

X9ZZ 

- 

- 


20 

0.7874 

6 

0.2362 

1590 

357 

46 000 

0.3 

0.012 

0.0076 

0.030 

AY 9 

- 

AY9ZZ 

- 

- 


24 

0.9449 

7 

0.2756 

3 710 
835 

36 000 

0.3 

0.012 

0.014 

0.031 

609 

609 Z 

609 2Z 

609 RSI 

609 2 RSI 


26 

1.0236 

8 

0.3150 

4 620 

1 040 

32 000 

0.3 

0.012 

0.020 

0.044 

629 

629 Z 

629 2Z 

629 RSI 

629 2 RSI 

10 

0.3937 

19 

0.7480 

5 

0.1969 

1 170 
263 

48 000 

0.3 

0.012 

0.0054 

0.012 

X10 

- 

X10ZZ 

- 

- 


#> 


19 

0.7480 

5 

0.1969 

1 480 

333 

38 000 

0.3 

0.012 

0.0055 

0.012 

61800 

- 

61800 2Z 

- 

- 

22 

0.0661 

6 

0.2362 

1 900 
427 

38 000 

0.3 

0.012 

0.016 

0.035 

61900 

— 

61900 2RZ* 

- 

61900 2RS1 

26 

1.0236 

e 

0 3150 

4 620 
1 040 

34 000 

0.3 

0.012 

0.0055 

0.012 

6000 

6000 Z 

6000 2Z 

6000 RSI 

6000 2RS1 

30 

1.1811 

9 

0.3543 

5 070 
1 140 

30 000 

0.6 

0.024 

0.032 

0.071 

6200 

6200 Z 

6200 2Z 

6200 RSI 

6200 2 RSI 

35 

1.3780 

11 

0.4331 

8 060 
1 010 

26 000 

0.6 

0.024 

0.053 

0.11 

6300 

6300 Z 

6300 2Z 

6300 RSI 

6300 2RS1 


t) 


1 This relers lo oil lubrication and moderate load. Consult SKF for lower ratings applicable to grease lubrication. 
Series 16100 — 16101, 16002 - 16072, also available. 

Series 6200 through 6220 and 6303 through 6317 are also available as precision bearings {ABEC 5). 

* Seal and shield versions 1mm wider than listed. 

3 Suffix 2RZ denotes rubberised shield. ^ 
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Material Selection 

The driving properties which pertain -to the Darn'ew 
blades are found below in the criteria. 

Preference Weightings 

0 - No difference. * 

1 - Very slightly more 
important-. 

2 - Slightly more 
important. 

3- "Reasonably more 
important. 

A l 1 I r _ 


ffllTFRlA 

8 

c 

b 



Tensi le to A 


h/s 

A 4 


A- 1 

E B 

Vi 

m 


TO" 

>. c 

■M 


£-3 

Fatigue d 


D = 3 



M- + 


r 


Materia l Selection 

The drjvina properties which pertain -fa -Hie 
Shaft are founa below 1Y1 4 he criteria. 


irrif>or+ftrVL. 

A- Much more 
imporUnf. 
5- Extremely more 

M I«sis MATRIX '"*•*"*■ 


u 


Toron A I umi num is -Hoe best- material far 
4he Shift. 



Phfferetr e Weigh!) pas 
0- liodiUenmt J 


1- Verysliohfly more 

impor+aht. 

2- Sliah+ly more 
- J i mpor+anf. 

b~ KW&nablv more 



PA-S 


Materia l Selection 

Thednvim properties which pertain to -the 
tpromill blades arefuind below in -the enter id. 

tVeference Weighty 

/’mi-,-,-,, » „ 0- No difference. 

I - Very slightly r» 


MV 

A- 10 

0 

3 

Sf V 

r C- 1 

b 0 

b-- o 

E 

£- 1 


3 ' KGciS 6 >nably rrior 
innpor+aru. 
•4 - M uch more 


-Analysis Matrix 


Steel 
A I uminuiTi 
T Itanium 
Magnesium 

boron Aluminum 
(evlar Woven 


/rt/Q. 
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Mb 

Vs 
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jL n 

V 2 

Mo 
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V 20 
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Lall-Xy - 400MPU 
VS. 


Yield : 025 MP« 
-Shear; 400MF& 

• 4 » 29/ COl 5 
MPa = 80 MFb 




/W r -%— = !*! & - = 9.C6N) ■. |.I3 xIO""Vm* 

” flOMFn goMFu eoMfti 
= I.ISxio^m^ /ioafti my = .00113mm 2 
U m ) 

Are«. = TLd 2 - .021 10mm 1 
^ d 2 * 1. 44 x lO-^mM 2 
d - 0.058 mm 

d = QJ)5 mm used |or | } • and rivet . 

check bendma mcmenf CT^ = ' ; 

; M=±uJl a »_L C2K/JC.oo2i)^<<(.4Uin' 

T 22 t m 

c^.oo 0025m 

1= 1 mR z = 1 (1. 0X I o^X 009025) 2 - 5 6>25X)0~ 17 
T= \9(pOOCO = 2 MPa 

cannot exceed 825 MP^ 

Weiahfe-- 


pin will hold 


im 

~Q Volume-- Trr’L = TT(.OO26cm) l (.03cm j )=l.O3x| O - 5 cm 5 
m - - p-V * 4.^/ (I63xi6-gcrw & )-- 7.2 x|q-% 

Volume: je i w H - fccmXacMX^cm)- - 4 . 0 cm J 

W : * ^ (.4.& crY |i) = 21. 22 c j. 

Tohxl : 21.220072^. X 4 pieces = 84, 6gP288 » k | 


i 
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T„ir d OOMPa - 0 OWPa 

^ - V 2 = 56^2. - ,^^7xiyV/ lOOm h^. 0035mM z 

"ft 80MPO 0DMftt \ IM J 

A™ - IL d ' 1 = . do^mm* 

d 1 -- 4.418*1 O' 3 mm* 

d - 0.0 b 1 mm 

d - 0. 61 -T- — /h3 our r^r 

s/tc or /? r A} 

Wel ^ Volume 4 pin: irr 2 L - TY(.ai 35 cm) 4 j. 5 em) 

r = 5.17 X I0' 5 cm 3, .4 

wf - p • V = 4 . 42 q/ ( 5 .~nxi(T 5 cm*r 2.55x10 . 3 
\ J/ m b 

To-hd •• , 

4 pins @ 2.55x10-^ » I 02xio^ 
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